Dispersion, controlled dispersion, and 
three applications 



o 

(N 



(N 



CZ2 

O 

-I— > 

q 

O 

• *— i 

Oh 



> 

in 

in 

o 
o 






by 



Douglas H. Bradshaw 

B.S., Applied Physics, Brigham Young University, 2000 
M.A., Biophysics, Johns Hopkins University, 2002 



DISSERTATION 

Submitted in Partial Fulfillment of the 
Requirements for the Degree of 

Doctorate of Philosophy 
Optical Science and Engineering 

The University of New Mexico 
Albuquerque, New Mexico 



May, 2010 



©2010, Douglas H. Bradshaw 



in 



Dedication 



To Gina 



IV 



Acknowledgments 



I have had the good fortune of learning from and studying under four exceptional 
mentors. Under Rudolph Wolfgang, I learned about the nature of optical pulses, 
about spectrometers and the power of seeing the world in reciprocal space. Under 
Ivan Deutsch, I learned some of the fundamentals of quantum mechanics and quan- 
tum optics. Under Michael Di Rosa, I benefited from a rare combination of freedom 
and involvement, along with encouragement to take my ideas and run and make my 
own mistakes. Finally, from Peter Milonni I have been slowly learning what it means 
to be a theoretical physicist who focuses on optics. 

Much of the most interesting physics for me happens not at a computer or in 
the lab but with a colleague near a chalk board. I'd like to acknowledge many 
pleasurable and fruitful discussions with Denis Seletskiy at the University of New 
Mexico. Our talks in the "Center for Amature Studies" form great memories around 
several interesting ideas. I'd like to acknowledge many similarly engaging discussions 
at Los Alamos National Laboratory with Michael Di Rosa, John Mcguire, and Peter 
Milonni. 

I've already mentioned Gina in my dedication, but she is the part of this whole 
experience that keeps on amazing me when I think about it. Somehow, Gina, you 
had the courage to embark on this long journey with me and our little growing family. 
I don't think that either of us realized going into this how long a Ph.D. can actually 
take. But your patience, your faith, and your persistant encouragement continue to 
amaze me. 

Finally, much of this dissertation was written while I was supported by a Los 
Alamos Laboratory Directed Research and Development grant. I am genuinely grate- 
ful for this support. 



Dispersion, controlled dispersion, and 
three applications 



by 
Douglas H. Bradshaw 

ABSTRACT OF DISSERTATION 

Submitted in Partial Fulfillment of the 
Requirements for the Degree of 

Doctorate of Philosophy 
Optical Science and Engineering 

The University of New Mexico 
Albuquerque, New Mexico 

May, 2010 



Dispersion, controlled dispersion, and 
three applications 

by 
Douglas H. Bradshaw 

B.S., Applied Physics, Brigham Young University, 2000 
M.A., Biophysics, Johns Hopkins University, 2002 

Ph.D., Optical Science and Engineering, University of New Mexico, 

2010 



Abstract 

Causality dictates that all physical media must be dispersive. (We will call a medium 
dispersive if its refractive index varies with frequency.) Ordinarily, strong dispersion 
is accompanied either by strong absorption or strong gain. However, over the past 
15 years several groups have demonstrated that it is possible to have media that 
are both strongly dispersive and roughly transparent for some finite bandwidth. In 
these media, group and phase velocities may differ from each other by many orders of 
magnitude and even by sign. Relationships and intuitive models that are satisfactory 
when it is reasonable to neglect dispersion may then fail dramatically. 

In this dissertation we analyze three such cases of failure. Before looking at the 
specific cases, we review some basic ideas relating to dispersion. We review some of 
the geometric meanings of group velocity, touch on the relationship between group 
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velocity and causality, and give some examples of techniques by which the group 
velocity may be manipulated. We describe the interplay between group velocity and 
energy density for non-absorbing dispersive media. We discuss the ideas of temporary 
absorption and emission as dictated by an instantaneous spectrum. We then apply 
these concepts in three specific areas. 

First, non-dispersive formulations for the momentum of light in a medium must 
be adjusted to account for dispersion. For over 100 years, there has been a gradual 
discussion of the proper form for the per-photon momentum. Two forms, each of 
which has experimental relevance in a 'dispersionless' medium, are the Abraham 
momentum, and the Minkowski momentum. If u is the angular frequency, n is the 
refractive index, h is Planck's constant, and c is the speed of light, then these reduce in 
a dispersionless medium to per-photon momenta of Hu>/(nc), and nhu/c respectively. 
A simple generalization of the two momenta to dispersive media entails multiplying 
each per-photon momentum by n/n g , where n g is the group refractive index. The 
resulting forms are experimentally relevant for the case of the Abraham momentum, 
but not for the Minkowski momentum. We show how dispersion modulates the 
displacement of a sphere embedded in a dispersive medium by a pulse. 

Second, pulse transformation in a nonstationary medium is modulated by the 
presence of dispersion. Dispersion may enhance or mitigate the frequency response 
of a pulse to a changing refractive index, and if dispersion changes with time, the 
pulse bandwidth must change in a compensatory fashion. We introduce an explicit 
description of the kinetics of dispersive nonstationary inhomogeneous media. Using 
this description, we show how the group velocity can modulate the frequency re- 
sponse to a change in the refractive index and how Doppler shifts may become large 
in a dispersive medium as the velocity of the Doppler shifting surface approaches 
the group velocity. We explain a simple way to use existing technology to either 
compress or decompress a given pulse, changing its bandwidth and spatial extent 
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by several orders of magnitude while otherwise preserving its envelope shape. We 
then introduce a dynamic descriptions of two simple media-one dispersive and one 
nondispersive. We compare the transformation of basic quantities like photon num- 
ber, momentum density, and frequency by a temporal change in the refractive index 
in a specific non-dispersive medium to those wrought by a temporal change in the 
group refractive index in a specific dispersive medium. The differences between to 
media are fundamental and emphasize the salience of dispersion in the study of 
nonstationary media. 

Finally, we note that the nature of a single optical cavity quasimode depends 
on intracavity dispersion. We show that the quantum field noise associated with a 
single cavity mode may be modulated by dispersion. For a well-chosen mode in a 
high-Q cavity, this can amount to either an increase or a decrease in total vacuum 
field energy by several orders of magnitude. We focus on the "white light cavity," 
showing that the quantum noise of an ideal white light cavity diverges as the cavity 
finesse improves. 
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Chapter 1 



Introduction 



The concept of group velocity has an interesting role in physics. On the one hand, 
its importance has been appreciated for a well over a century. For example, group 
velocity was lucidly described by Lord Rayleigh as early as the late 1800s [1, 2]. In 
a letter to Nature entitled "On the velocity of light" he explained [2], 



It is evident however that in the case of light, or even of sound, we 
have no means of identifying a particular wave so as to determine its 
rate of progress. What we really do in most cases is to impress some 
peculiarity, it may be of intensity, or of wave-length, or of polarization, 
upon a part of an otherwise continuous train of waves, and determine 
the velocity at which this peculiarity travels. Thus in the experiments of 
Fizeau and Cornu, as well as in those of Young and Forbes, the light is 
rendered intermittent by the action of a toothed wheel; and the result is 
the velocity of the group of waves, and not necessarily the velocity of an 
individual wave. 
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Group velocity also played a critical role in de Broglie's theory of matter waves [3]. 
While the phase velocities of de Broglie waves associated with a moving object are 
superluminal, de Broglie showed that the group velocity of the waves matches the 
velocity of the object, relating the group velocity to the propagation of mass/energy 
through space. The important role that group velocity and group index play in the 
performance of spatial interferometers was also recognized generations ago [ ]. 

On the other hand, differences between the group velocity and the phase velocity 
are often neglected. In optics, this is often reasonably practical. Strong dispersion is 
generally accompanied by strong absorption or strong gain, while transparent media 
tend to be only weakly dispersive. Since the bulk of optical experiment and theory 
(sensibly) involves light traveling through roughly transparent media, dispersion is 
often negligible. We can illustrate this with a numerical example. Suppose that 
we are concerned with an etalon or an optical cavity that consists of BK7. At a 
vacuum wavelength of 633 nm, BK7 has a refractive index of 1.515 and a group 
refractive index of 1.537. The formula for the free spectral range represented in some 
standard introductory optics texts (see, for example, Verdeyen or Guenther [5, 6]) 
is Ai>fsr — c/(2nd), where n is the refractive index and 2d is the round-trip cavity 
length. A formula that takes first order dispersion into account (but in its turn 
neglects higher orders of dispersion) would be Aufsr = c/(2n g d), where n g is the 
group index. For BK7 at 633 nm, the expressions differ from each other by 1.4%. 
From this perspective, several results of the last decade or so are interesting because 
they combine strong dispersion with reasonably small absorption for select, narrow 
bandwidths. In 1999, Hau, Harris, Dutton, and Behroozi reported that they had 
succeeded in slowing the propagation of a pulse of light to a mere 17 m/s [ ]. If this 
medium could be placed in a cavity then neglect of dispersion could lead to a larger 
error in the estimate of the free spectral range than in the case of BK7-instead of 
an error of 1.4% it would lead to an estimate that is too large by a factor of over 17 
million for the relevant wavelengths. 
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The 1999 demonstration of slow light in a weakly absorbing media was just one in 
a string of experiments showing increasing control over dispersion. In 2000, Wang, 
Kuzmich, and Dogariu demonstrated nearly lossless pulse propagation through a 
medium with negative group velocity, implying that the peak of a pulse left the far 
side of the negative group velocity medium before the incoming pulse peak actually 
entered the medium [8]. This corresponds to a negative group index. Finally, in 2001 
Liu, Dutton, Behroozi, and Hau reported that they had completely stopped a pulse 
of light and then had recovered it and sent it on its way [ ]. This corresponds to 
"stopped" light that is later "started" again. It not only demonstrated that the group 
velocity could be brought to 0, but that it could be controlled dynamically while a 
pulse was propagating through the medium. All three of these feats, those of slow, 
fast, and stopped light, were performed in the absence of substantial absorption or 
gain or pulse deformation. Together, they represent newly opened physical regimes 
where dispersive effects in roughly transparent media are not only non-negligible, 
but may dominate the answers to particular questions. The opportunity represented 
by the opening up of these new regimes is a primary motivation for this dissertation. 



1.1 Dissertation outline 

In this dissertation we seek to do two things. First, we lay a foundation for un- 
derstanding the physics of dispersion. In Chapter 2, we review the definition of 
group velocity. We look at some of the geometric meanings of group velocity. We 
distinguish the group velocity from the phase velocity, the signal velocity, and the 
energy velocity. We also briefly review some of the work that has been done on the 
relationship between anomalous dispersion and causality. Following this we give a 
brief review of some of the techniques whereby strong dispersion with minimal gain 
and/or absorption is brought about. In Chapter 3, we look at the interplay between 
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energy density and dispersion. Following the work of Peatross et al. on the role 
of the instantaneous spectrum in pulse propagation through dispersive media [10], 
we also explore the relationship between dispersion and temporary absorption and 
emission in "absorptionless" media. 

Second, we apply these ideas to three specific questions. The relationship between 
dispersion and absorption requires that dispersion have an impact on momentum 
exchange between field and medium. In Chapter 4, we explore the relationship 
between dispersion and photon momentum, specifying how the nondispersive forms 
for the Abraham and Minkowski momentum densities are altered when dispersion 
is taken into account. We predict that by modulating the group velocity of a pulse 
in a medium it should be possible to modulate the displacement of a small particle 
embedded in that medium. 

One interesting aspect of controllably dispersive media is that in many cases the 
optical properties of a medium may be altered while a wave is propagating through 
it. This suggests that it may be profitable to look at controllably dispersive media 
from the perspective of nonstationary electromagnetics. In Chapter 5, we ask what 
effects dispersion can have on the transformations wrought upon a pulse traveling 
through a nonstationary medium. We begin with a general kinetic approach based 
on the preservation of discrete translational symmetries. Using this approach we find 
that dispersion may modulate the frequency response of a pulse to temporal changes 
in refractive index, and derive simple analytical expressions for Doppler shifts in 
dispersive media when group velocity dispersion may be neglected. These show 
that dispersion may lead to large reflective Doppler shifts when the group velocity 
grows close to the velocity of a moving surface. Using the kinetic formalism, we also 
show how temporal and spatial control of the group velocity of a medium may be 
combined to scale pulse duration (and so pulse bandwidth and longitudinal extent) 
over many orders of magnitude without otherwise altering the pulse envelope. This 
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explains and generalizes a recent numerical proposal for the compression of pulses 
using magnetized dispersive nonstationary plasmas [] L]. We then supplement these 
general considerations with a comparison of two specific transformations, one by a 
temporal change in the refractive index in a non-dispersive medium, and the other 
by a temporal change in the group index in a dispersive medium. In order to make 
a full comparison, we introduce boundary conditions for each transformation. For 
the non-dispersive case, the boundary conditions are taken from a classic work by 
Morgenthaler [12]. For the dispersive case, we use an idealized version of the bound- 
ary conditions that were apparent in the stopped-light experiment of Liu, Dutton, 
Behroozi, and Hau ['.)]. We compare the effect of these two different transformations 
on 20 different fundamental quantities associated with a propagating pulse. 

In Chapter 6, we address the effect of anomalous intracavity dispersion on the 
quantum field noise associated with a single (pseudo-)mode in an optical cavity. 
We do so using two approaches. The first approach follows earlier work by Drum- 
mond [13] and Milonni [J l]. In this approach we assume a lossless dispersive cavity 
and show via the classical dispersive energy density that the electromagnetic field 
strengths all scale with the square root of the group velocity. Associating a harmonic 
oscillator which each mode leads to noise terms whose associated fields must similarly 
scale with the square root of the group velocity. When the dispersion is anomalous, 
the classical expression for the energy density implies that the energy associated with 
the fields alone is larger than the total modal energy. This is because the total modal 
energy also includes a term that represents temporary absorption by the medium. 
In lossless, anomalously dispersive media this term is negative, implying temporary 
gain, in accordance with the dictates of the instantaneous spectrum. In practical 
terms, this means that the medium energy is negative because it has temporarily 
donated energy to the electric field. When the average group index approaches zero 
as in the "white light cavity" ], the total classical energy density goes to zero 
unless the field energy goes to infinity. Thus the field energy associated with even 
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the vacuum state of such a cavity diverges. These considerations motivate us to 
calculate the quantum field noise of anomalously dispersive cavities, and particularly 
of white light cavities, from a second perspective. We use a simple approach where 
the pseudo-modes of an open cavity may be probed by true modes or "modes-of- 
the-universe [16]." This second approach allows us to corroborate the results of the 
first without relying on an expression for the electromagnetic energy density of an 
anomalously dispersive medium. Because this approach allows for an open cavity, it 
not only reveals the amount of vacuum field noise associated with a particular reso- 
nance but also shows the effect of dispersion on spectral distribution of that noise. 
We first show analytically that the second approach leads to the same conclusion as 
the first approach if we assume perfect finesse. We then relax this assumption and 
show numerically in a specific physical model how the quantum field noise depends 
on finesse. The divergence associated with a white light cavity for the case of infi- 
nite finesse is relaxed when reasonable physical assumptions are made. However, the 
quantum field noise associated with an anomalously dispersive cavity mode remains 
substantially larger than that of the corresponding evacuated cavity. 

These last three paper-like chapters represent three papers. Chapter 4, on mo- 
mentum, is taken from a paper [17] that I coauthored with Zhimin Shi and Robert 
Boyd (University of Rochester), and with Peter Milonni (Los Alamos and University 
of Rochester), who was the primary author. Professor Boyd's group is currently 
planning experiments to test some aspects of the theory developed in that paper, 
which has been published in a special issue of Optics Communications in memory of 
Krzysztof Wodkiewicz. The work in this paper has been cited in two manuscripts by 
Rodney Loudon, FRS, and Stephen Barnett, FRS, to be submitted for publication 
in Physical Review Letters and Proceedings of the Royal Society of London. Chapters 
5, on pulse transformations by a dispersive nonstationary medium, and 6, on the in- 
terplay between dispersion and the vacuum field energy associated with a particular 
cavity mode, are papers which I have coauthored with Michael D. Di Rosa. I am the 
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Journal Reference 


Coauthors 
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Optics Communications Volume 283, Pages 650-656 
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Ready for submission 
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Ready for submission 
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primary author of these last two chapters. 



Chapter 2 

Group velocity, the group 
refractive index, and controlled 
dispersion 



When we use the word "dispersive" to describe a medium, we mean that the refractive 
index of that medium varies with frequency. Another way to say this is that in a 
dispersive medium the group velocity differs from the phase velocity and the group 
index differs from the refractive index. In order to understand the implications of 
controlled dispersion, we need to understand the meaning of the group velocity and 
of its companion, the group index. 

In a letter to Nature entitled "On the velocity of light," Lord Rayleigh wrote [2], 

I have investigated the general relation between the group velocity U and 
the phase velocity V. It appears that if k be inversely proportional to 
the wavelength, 

rr _ dim 

dk ' 
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and is identical with V only when V is independent of k, as has hitherto 
been supposed to be the case for light in vacuum. 



The definition of the group velocity employed by Lord Rayleigh in 1877 and 1881 
is the one we use today. Writing k = nu/c and V = c/n, we write the group velocity 
in its more standard form, 

du 



The group index is defined analogously to the phase index: 



n g = -= c— . 2.2 

v n du 



The purpose of this chapter is to explore some of the fundamental ideas associated 
with the group velocity, the group index, and their control. 



2.1 Some geometric meanings of group velocity 
and the group index 

The definition of the group velocity as given in Eq. (2.1) is really a geometrical 
definition. It says that the group velocity is the ratio of infinitessimal changes in 
temporal periodicity to infinitessimal changes in spatial periodicity. In this section 
we review some simple geometrical consequences of dispersion that have practical 
optical consequences. 
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2.1.1 The group velocity is the velocity of fixed phase re- 
lationships between copropagating plane waves with 
infinitesimal spectral separation 

Electromagnetic disturbances propagating through a lossless, homogeneous, isotropic, 
linear medium may be decomposed as a summation of plane waves. The behavior of 
the disturbance over time depends on the interference pattern formed by the waves, 
which is dictated by the relative phases between the many components at each point 
in space and time. If two plane waves have phases 0i = fax — wi£ + #i and 02 = fax — 
oj2t+02, the relative phase maybe taken to be 02— 0i = (fa — ki)x— (o>2 — u)i)t+02— 9i- 
The velocity (v 5) of a beat between these two waves is simply the velocity of a plane 
of constant relative phase. This velocity is given by 

Vb = 1 — r- ( 2 - 3 ) 

fa - fa 

Taking the limit for the beat velocity as the waves grow closer to each other in 
temporal and spatial frequency yields the group velocity. To the extent that the beat 
velocity is constant for the phase relationships between all of the salient frequency 
components of a pulse, the pulse will propagate without distortion. 

The velocity of fixed phase relationships between plane waves with finite 
spectral separation may be found through an averaged group index 

Interestingly, for finite frequency differences, the beat velocity is not that given by 
the spectral average of the group velocity between the two frequencies. Rather, it is 
determined by the a spectral average over the group index, where the group index is 
defined as 

% = c-. (2.4) 
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We can see this by writing h 2 in terms of h\ as 



h = h + I —du = hi + / n9 ^' duo. 
Im duJ J ui c 



Thus, 



k 2 - hi = (u> 2 -u>i) — , (2.5) 

c 

where the line over n g is used to denote a spectral average, so that 



f^ 2 n a (u)du 
u 2 - Ui 

Thus, the beat velocity may be written in terms of the spectrally averaged group 
index as 

k 2 - h n 9 ' 
This definition, based on a spectral average over all frequencies between the frequen- 
cies of interest, accounts for group velocity dispersion as well as all higher orders of 
dispersion. If these terms are negligible for the bandwidth of interest, n g may be 
replaced with n g . 

A pulse will propagate through a medium without distortion only to the extent 
that Vb is similar for each pair of frequency components. 



Relating the velocity of fixed phase relationships to the velocity of pulse 
propagation in a lossless medium. 

Having established that n g and v g can be interpreted in terms of the speed of constant 
phase relationships, we now explore the relationship between the speed of a fixed 
phase relationship and the speed with which a pulse will travel through a medium. 



11 



Chapter 2. Group velocity, the group refractive index, and controlled dispersion 

Briefly, the velocity of a pulse is generally considered to be the velocity of its 
envelope. The envelope of a pulse is defined by the phase relationships between the 
frequency components of that pulse. Thus, the velocity of those phase relationships 
is the velocity of the pulse. 

We now make this argument in more detail. First, we introduce the group and 
phase velocities as limits of two velocities that naturally arise when different fre- 
quency components are summed. In doing so, we find that although these velocities 
are, in practice, applied to pulses with finite bandwidths, they acquire exact mean- 
ings only over spectra sufficiently narrow that changes in n (for the case of the phase 
velocity) and n 9 (for the case of the group velocity) can be neglected across the 
relevant range. We then find that if these changes can be neglected then the pulse 
envelope is defined by the phase relationships of the frequency components of the 
pulse. Thus, the speed of these phase relationships defines the speed of the pulse. 

The phase and group velocities are only well-defined over frequencies 
ranges sufficiently narrow such that n and n g can be considered to be 
constant. A pulse can be thought of as a sum of plane waves. If we consider prop- 
agation in a single dimension, we can go further and consider a pulse in terms of a 
sum of frequency components. To begin to see how the group velocity works, it is 
sufficient to start with the very simple example of a sum of two sinusoids of equal 
amplitude. We might represent such a sum as 

E = cos(cji£ — k\x + cf>i) + cos(w2t — k^x + $2). (2-8) 

We can rewrite this sum as a multiple of two sinusoidal terms by using the identity 

fA + B\ (A-B 
cos(A) + cos(.B) = 2 cos I — - — cos 



This gives 

E = 2 cos (out -kx + (p) cos (Auut - Akx + A0) , (2.9) 



12 



Chapter 2. Group velocity, the group refractive index, and controlled dispersion 

where u = (u 2 + Wi) /2 and Au = (o>2 — U)i) /2 and similarly for the fc and terms. 

If CU2 and u/i are similar in size, to each other, the difference between them (2 Au) 
is much smaller than their sum (2a)). The slowly varying term becomes an envelope 
gradually damping out and reviving the oscillation of the second term. 

The velocity associated with a constant value for the rapidly oscillating term is 
iu/k, while that associated with the slowly oscillating term is Au/Ak. Rewriting 
these terms by expanding them in terms of values for u and n gives 

v f = -. (2.10) 

UJ2U2 C + UiTli C 



and 

u 2 - Ui 



(2.11) 



u 2 n 2 /c + uifii/c 

where Vf is the velocity of the fast oscillations and v s is the velocity of the slow 
oscillations. While these expressions are related to the group velocity and the phase 
velocity, they are not exactly equal to them. To relate them to the group and phase 
velocities, we must either make assumptions about n(u) or about the frequencies. 
Taking the limit of the two velocities as (u 2 — cui) — > gives Vf = c/n = v p and 
v s = c/n 9 = v 9 . 

Here we see clearly that the group and the phase velocity are only well defined as 
functions of frequency. Since any pulse must have components taken from a range of 
frequencies, these velocities as applied to that pulse are approximate to the extent 
that n and n g change across those frequencies. 

Changes in the envelope of spectrally narrow pulses in lossless media 
can come only from changes in phase relationships between its frequency 
components. A pulse can be thought of as a sum of sinusoids. That is, 

E(t) = f duE(cu)e- iwt dco, (2.12) 
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where E{uS) = \E{u)\e~ l ^ u \ If the amplitudes are constant then only the phases 
change. The absolute phase determines the position of the envelope as a whole. The 
relative phases determine the shape of the envelope. 

The group index and the free spectral range of an optical cavity 

In Chapter 1 we used the free spectral range of an optical cavity as an example 
of a case where dispersion is generally ignored but can have real consequences. A 
commonly used form [5, 6] for the free spectral range is 

Q 

Au F SR non = ——, (2.13) 

2Ln 
where 1L is the round trip cavity length. We asserted that a more correct form for 
the free spectral range would be 

AvFSRdis = wj— . ( 2 - 14 ) 

but acknowledged that this is also an approximation that ignores higher order dis- 
persion. 

However, we are now in a position to argue that the difference between Eqs. (2.13) 
and (2.14) is more fundamental in nature than the difference between Eq. 2.14 and 
forms that account for higher order dispersion. 

The round trip phases associated with two neighboring cavity resonances must 
differ from each other by 2tt. Thus, 

(k n+1 -k n )(2L) = 27r. (2.15) 

We can relate the difference in wave numbers in Eq. (2.15) to a difference in frequency 
using Eq. (2.5), giving 

2ttc 
2n g L 
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so that the free spectral range, Aufsr, is given by 

&VFSR = ^7, (2-16) 

2n g L 

where n g is given by Eq. (2.6). This expression is exact. Thus, accounting for higher 
orders of dispersion leads to taking a spectral average of the group index between the 
two neighboring resonances. On the other hand, n g has no fixed relationship with 
the spectral average of the refractive index between the two resonances. 

In summary, there is a difference between the kind of approximation we make 
when we ignore dispersion and the kind of approximation we make when we ignore 
only group velocity dispersion and higher orders of dispersion. If we ignore dispersion, 
the formula for the free spectral range becomes Aufsr ~ c/2nd. Taking only first 
order dispersion into account, we get Aufsr ~ c/2n g d. Finally, we have just seen 
that if we account for all orders of dispersion we get Avfsr — c/2n g d. As a cavity 
grows longer, the free spectral range becomes smaller and n g — > n g whether or not 
higher orders of dispersion are significant. In a dispersive medium, n g will not make 
a similar approach to n regardless of cavity length. Intuitively, we can understand 
the more fundamental relationship between n g and the spectral separation between 
neighboring resonances in terms of the defintion of the group index, n g = c(dk/du). 

2.1.2 Group optical path length 

A pulse that propagates through a spatial interface from a medium with one group 
velocity to a medium with a second group velocity is scaled longitudinally [ ] but 
not temporally (when a pulse crosses an interface with a temporal component, for 
example a moving interface, this is no longer true-see Chapter 5). If the pulse dura- 
tion is given by r p and the longitudinal extent of the pulse by I then (in the absence 
of distortion) the two quantities may be related by the group velocity according to 
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r p = l/vg = rigl/c. Just as the pulse duration is conserved at a stationary interface, 
so is the product n g l. 

Following Candler [ i], we now extend the idea of the optical path length to define 
a group optical path length. The optical path length associated with a distance d in 
a homogeneous medium is just the distance multiplied by the appropriate refractive 
index, 



J OP 



nd. (2.17) 



The phase accrued by a wave in traveling a distance d in a medium with a refractive 
index n is the same as that accrued by a wave of the same frequency traveling through 
a vacuum distance nd. 

The group optical path length associated with a distance d in a homogeneous 
medium is just that distance multiplied by the appropriate group index, 

Loop = n g d. (2.18) 

The difference between phases accrued by two infinitesimally different frequencies as 
they travel along a path with a group optical path length of Loop is the same as the 
difference they would accrue if they traveled an actual distance of Loop through the 
vacuum. In this sense a milimeter of propagation through a slow light medium may 
be equivalent to a kilometer of propagation through vacuum if n g = 10 6 . 

2.1.3 Group optical path length and spatial interferometers 
used as spectral filters 

Many interferometers that are used for spectral separation actually depend on wave- 
length. Fabry-Perot interferometers, Young's double slit experiment, and diffractive 
gratings, for example, may all be analyzed in terms of wavelength. If the mapping 
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between wavelength and frequency is altered, the frequency response of such inter- 
ferometers is altered while the wavelength response is unchanged. To emphasize this 
distiction, we will refer to such interferometers as spatial interferometers. 

The interference patterns of spatial interferometers depend directly upon wave- 
length and only indirectly upon frequency. The spectral resolving power of a spatial 
interferometer comes from phase differences in optical path lengths for different col- 
ors. That difference is related to the group optical path length: 

d{kL) dk rig 

duo duj c 

Thus, as was pointed out by Candler [ 1], the spectral resolving power of spatial filters 
is conveniently given in terms of the group optical path length. We apply this idea 
to an optical cavity in greater detail in Chapter 6 and particularly in Appendix B. 

A second way to understand the function of the group index in spatial inter- 
ferometers is by thinking of the group index as a map between relative changes in 
wavelength and changes in frequency. To see this more clearly, we can rewrite the 
group index as 



d In A , 

n g = -n- 2.19 



(see Appendix A). In vacuum, n g — 1, and n = 1, meaning that a .01% increase 
in frequency leads to roughly a .01% decrease in wavelength. In a moderately slow 
medium, n g might be 100 while n remains at approximately 1. In this case a .01% 
increase in frequency leads to roughly a 1% decrease in wavelength. Thus, with no 
change to the spatial resolution of an interferometer, its spectral resolution may be 
enhanced 100-fold. 
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2.2 Group velocity, information velocity, energy 
velocity, and causality 

Lord Rayleigh argued that in measurements of the speed of light, it is not typically 
the phase velocity that is measured, but the velocity of a peculiarity or a change 
impressed on that light [1, 2]. That peculiarity would travel at the group velocity 
Later, in de Broglie's thesis on the wave nature of matter, de Broglie derived super- 
luminal phase velocities for the component waves of a moving particle. Referring to 
Rayleigh's work, he showed that although the phase velocities associated with his 
matter waves would be superluminal, the group velocity of those waves would match 
the velocity of the particle [3] . 

From Lord Rayleigh's work comes the idea that the group velocity may be used 
as a signal velocity. From de Broglie's work we see that the group velocity can 
function as a velocity of energy transport. However, equating the group velocity 
with either the signal velocity or the energy velocity becomes problematic when the 
group velocity exceeds the speed of light. Superluminal group velocity occurs even in 
ordinary media; the group velocity at the spectral center of an absorbing resonance 
is generally superluminal (see Section 2.3). 

The subtleties associated with superluminal group velocity have been grounds for 
careful thought for over a century and still drive a large portion of the discussion of 
dispersive media; a valuable dissertation (or perhaps several) could still be written 
on the topic. However, because the relationship between causality and dispersion is 
beyond the scope of this work, we make just a few points and refer the reader to [19 
for a more detailed introduction. 

First, causality is built into the response functions of the media used to obtain 
superluminal group velocities. Thus, the question is not whether causality is broken 
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by such media-it cannot be. Rather, the goal is to understand how causality is 
upheld. 

Second, the task of understanding how causality is upheld is not as simple as it 
may first seem because it requires that we develop precise definitions of quantities 
like the information velocity and the energy velocity. In many practical situations 
the group velocity works perfectly well as a proxy for each of these. However, the re- 
ality of absorptionless anomalous dispersion shows that both velocities are in general 
independent of the group velocity. 

Third, the difficulty of the problem is in fact an opportunity to come to a bet- 
ter understanding of what we mean when we reference the energy velocity or the 
information velocity of a disturbance. This is the lens through which the rest of the 
points in this section should be seen. 

2.2.1 Information velocity 

Around 1910, the question of causality in the presence of superluminal group velocity 
was the subject of several conferences. Responding to the topic, 

Sommerfeld demonstrated theoretically that the velocity of the front of a 
square-shaped pulse propagating through any medium is identically equal 
to c and hence relativistic causality is preserved. In a follow-up study, 
Brillouin suggested that the group velocity is not physically meaningful 
when the dispersion is anomalous because the pulse becomes severely 
distorted [20]. 

An English summary of this part of the history can be found in Brillouin's Wave 
propagation and group velocity [21]. 
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Brillouin's suggestion that a pulse will become severely distorted as it propagates 
under anomalous dispersion is likely to be correct if the distance through which 
a pulse propagates under anomalous disperion is sufficiently long. However, it is 
demonstrably incorrect for propagation over shorter distances. Garrett and Mc- 
Cumber showed via simulations that a pulse passing through a finite anomalously 
dispersive medium may retain its overall shape while its average position advances 
faster than the speed of light for a limited distance [: ] . Later analysis and experi- 
ment have demonstrated superluminal pulse propagation without gross distortion in 
active media [23, 24, 25, 8, 20]. Pulse distortion does not therefore provide satisfac- 
tory explanation for the preservation of causality. 

Sommerfeld's statement on the propagation of the front of a square-shaped pulse 
leads us to consider what it means to have information communicated by the electro- 
magnetic field. In his simple model, he associates that information with the moment 
at which a pulse is turned on. The spectrum associated with this instant transition 
is infinite, and any physical response function must therefore dictate that the signal 
associated with this moment moves at c. 

This concept has been generalized by Garrison et al [ ], who propose that a 
signal is associated with a point of non-analyticity. Because any such point will have 
an infinite spectrum, it will experience a total group index of 1 and move at c. If the 
information velocity can be associated with such points, then our problem is solved. 

However, this definition of group velocity leaves us without a way to understand 
the superluminal transmission of a smooth pulse. It may be impossible to generate a 
pulse that is smooth for all time because such a pulse would have an infinite extent. 
However, it certainly seems possible to generate, for example, a Gaussian pulse that 
is smooth for that portion of the pulse which has a power over some minimum 
threshold. Because of the fast decay of the wings of a Gaussian, that threshold could 
be very low without making the duration of the analytic region unreasonably long. 
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For such a pulse we may need a different definition for information velocity. 

Kuzmich et al. j ] used signal-to-noise ratios to analyze the faster-than-light 
propagation associated with single photon pulses propagating through a non-lossy 
superluminal medium. They showed that the speed of propagation of the photon as 
given by its expected arrival time at a detector could actually exceed c. However, 
the signal velocity, defined operationally as using the signal to noise ratio, actually 
decreased because of the anomalously dispersive medium due to the effect of quantum 
fluctuations. 

This slower signal velocity, as judged by the signal to noise ratio, was later verified 
in practice (see, for example, Stenner et al. [20]). 



2.2.2 Energy velocity 

Just as causality prevents a signal from travel faster than the c, it is also prevents 
energy from traveling faster than the speed of light. However, as we have just 
mentioned, Kuzmich et al. showed that the speed of propagation of the photon 
as given by its expected arrival time at a detector could indeed exceed c [ ]. 

Evidently, the fact that energy may not propagate faster than c does not mean 
that the energy associated with a particular wave cannot travel faster than light if 
the wave is not a closed system. That is, if energy is added to the front of a moving 
pulse and taken from the back of the pulse then the electromagnetic field energy 
associated with that pulse may have an average velocity that is faster than the speed 
of light. This is allowed through energy exchange with the medium through which 
the light is traveling. 

An analogy is useful here. There is a limit to how fast an orange may move 
around the globe . . . perhaps it is currently about the speed of a fast jet. However, 



21 



Chapter 2. Group velocity, the group refractive index, and controlled dispersion 

if your brother is in London and you live in Albuquerque, it is still possible for you 
to give him an orange faster than this so long as there are oranges in London. You 
must simply call someone in London and ask them to give your brother an orange. 
You can promise to send another orange to repay them following the next flight out 
of Albuquerque. 

One interesting point with regard to this metaphor is that if you are to give your 
brother the orange early, you must still be able to send the message somehow. That 
message cannot move faster than the speed of light. In the case of a pulse propagating 
through a superluminal medium, the message corresponds to the leading edge of the 
pulse, which must be amplified. The pulse can never pass this leading edge. 

The metaphor also brings the question of how the exchange of energy between 
pulse and medium is mediated. We will address this question in Section 3.3. 



2.3 Controlled dispersion 

Controlling dispersion means controlling the slope of the real part of the refractive 
index. Because the real and imaginary parts of the refractive index are linked via 
the Kramers-Kronig relation, we can alter the real part by changing the imaginary 
part. We will consider three media from the perspective of their linear electric 
susceptibility, which we will call %. 

In a linear medium, the complex refractive index, n, is given by n = ^e r fi r = 
y/(l + x)(l + Xm)i where Xm is the magnetic linear susceptibility. For x << 1 
in a nonmagnetic medium, the real and imaginary parts of x are related to the 
real and imaginary parts of the complex refractive index, rj, via rj r ~ 1 + Xr/2, 
and rji = x«/2- In this case we may use the electric susceptibility to conveniently 
quantify the dispersion and absorption of different media. A dispersive but roughly 
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Figure 2.1: At zero detuning (A = 0), a simple 2-level absorptive resonance combines 
strong anomalous dispersion and strong absorption. 



transparent medium will have a high ratio u(dxr/du)/xi 



2.3.1 Two-level resonance 

Figure 2.1 shows the real and imaginary susceptibilities, Xr and x%i associated with a 
probe weakly interacting with a 2-level absorptive transition. (In a sufficiently dilute 
medium, \x\ « 1 and the real and imaginary susceptibilities are proportional to the 
refractive index and absorption coefficients.) In Figure 2.1, the peak of absorption 
corresponds to strong negative dispersion at A = 0. The relationship between the 
real and imaginary parts of analytic function, as given by the Kramers-Kronig rela- 
tionship is roughly similar to that of a derivative. Thus, the peak of Xi corresponds 
to the maximum slope of Xr, while the peaks of Xr correspond to the maximum slope 
of Xi- Using this simple heuristic, it is easy to guess at the consequence of placing 
two absorptive resonances side by side. 
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Figure 2.2: In between two absorptive resonances placed side- by-side, there is strong 
normal dispersion and group velocity dispersion is minimized. 



2.3.2 Two absorptive resonances side by side 

As is represented in Figure 2.2, two absorption lines may be combined to form a 
slow-light medium with only moderate absorption. If this profile is combined with 
a second medium that exhibits broad-band gain, the overall absorption may be mit- 
igated without severely dampening the dispersion. As symmetry between the two 
absorption lines is approached, all even derivatives of \ r go to zero. However, this 
arrangement is no better than the wing of a single resonance in terms of maximizing 
the slope of the real part of the susceptibility while minimizing the magnitude of the 
imaginary part. 

2.3.3 Strong, controllable normal dispersion through Elec- 
tromagnetically Induced Transparency 

A more powerful approach utilizes Electromagnetically Induced Transparency (EIT) 

28, 2 ]. In this approach, diagrammed in Figure 2.3, quantum interference between 

two absorptive pathways leads to the potential for stronger dispersion coupled with 
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Figure 2.3: EIT. Together, a coupling beam (fl c ) and an atom form an effective 
medium as seen by the probe beam (Q p ). Destructive quantum interference between 
two alternative absorption pathways allows for minimal absorption in combination 
with strong dispersion at A = 0. 



less absorption. The EIT configuration has the additional advantages of more perfect 
symmetry (leading to a complete disappearance of even derivatives of Xr at A = 0) 
and controllability. While an ideal EIT system would have no absorption, in practice, 
loss due to decay, inhomogeneous broadening, spatial inhomogeneity, and various 
other sources of decoherence is always present. EIT configurations have been used 
to achieve vary slow and even vanishing group velocities [7, 9, 30]. 



2.3.4 Controllable anomalous dispersion via a double Ra- 
man scheme 

Although an anomalously dispersive counterpart to EIT has not yet been devised, two 
gain lines may be combined, yielding controllable anomalous dispersion via the dou- 
ble Raman scheme represented in Figure 2.4. The two gain lines lead to anomalous 
dispersion coupled with moderate gain. These may be combined with a broad-band 
absorber to yield a medium with no net gain but with strong anomalous dispersion. 
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Figure 2.4: At zero detuning (A = 0), a simple 2-level absorptive resonance combines 
strong anomalous dispersion and strong absorption. 



This technique has been used to obtain strong anomalous dispersion with minimal 
gain [8, 20, 31]. 
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Chapter 3 

Energy density and temporary 
absorption in lossless dispersive 
media 



In 1999, Hau and Harris published a paper that concretely described some of the 
features of a pulse that is introduced into a slow-group- velocity medium [18]. One of 
the most striking features of this transition is the spatial compression of the pulse. 
When a pulse leaves a vacuum and enters a slow-light medium, it may be compressed 
longitudinally by many orders of magnitude. However, in the case they examined, 
that of a pulse entering an EIT medium (see Section 2.3), the electric permittivity 
and the magnetic permeability in the slow light medium are very close to their 
vacuum values for fields near EIT resonance. Therefore the values of the squared 
electric and magnetic fields remain continuous across the boundary even as the pulse 
is compressed. The total field energy associated with the pulse then scales with the 
group velocity. When that velocity is slowed, very little of the original pulse energy 
remains in the field. Rather, it is coherently stored in the slow-light medium. 
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The partitioning of energy between field and medium is, in general, a subtle 
subject. However, in the special case where absorption may be ignored it becomes 
simpler. In this chapter, we examine the effect of dispersion on the energy density 
of a dispersive, lossless medium. 

3.1 Energy density in a lossless, dispersive medium 

Jackson [ _] credits Brillouin [21] for the original formulation of the energy density of 
a of a spectrally narrow wave traveling through a dispersive medium with negligible 
loss. We write this expression, considering the case where the imaginary portions of 
e and /i can be ignored, as 

B =i^M(E.E> + ^( Wo )(H.H>, (3.1) 

where the angled brackets denote a cycle average. 

In an isotropic medium, D = eE and B = /iH, and we can rewrite this expression 



{ f dine \ (D-E) + / + d\nu.\ (B • H) 



as 



dliaoj J 2 \ dlnojj 2 

For a plane wave propagating through a lossless medium, |H| = ^/e//j|E|, and /i and 
e are real so that 

(D-E) = (B-H). 

Using this fact and the identity 

n a 1 cHne 1 din a 

n 2dmui 2 d In /x ' 

which is valid when the imaginary parts of e and /i are negligible (see the end of 

Appendix A and Eq. (A. 19)), we can write 

u 
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But the term in brackets is just the cycle average of the energy density of an elec- 
tromagnetic field in a linear, nondispersive medium. We may therefore write 

n 
u = -%, (3.2) 

n 

where uq is the cycle-averaged energy density in the case where there is no dispersion, 
written 

(D-E) (B-H) , , 

«o = ^"^ + ^~ r L - (3-3) 

Note that the ratio between the magnitude of the Poynting vector and the energy 
density for a traveling quasimonochromatic wave in a dispersive medium is therefore 

ISI n IE x HI c ,„ A . 

L-L = —J L = — . (3.4) 

U Tig UQ Tig 

Expression (3.2) for the energy density in a dispersive absorptionless medium is 
intuitive and matches our expectations based on the 1999 paper by Harris and Hau 
[18], and comparisons between u and u$ lead naturally to a partitioning between 
energy that is stored because of dispersion and the total energy. When n « 1, as 
was the case in their paper, the energy can be neatly partitioned into a fraction, 
I /rig, that is stored in the form of field energy and a fraction, (n g — l)/n g , that 
is stored by the medium. When n ^ 1, there seem to be two different forms of 
storage by the medium, one that is related to dispersion, and another that is related 
to the phase velocity. We can form a partition between a dispersively stored energy 
fraction (((n g /n) — l) / (n g /n)) and non-dispersively stored energy fraction (l/(n g /n)). 
However, the non-dispersively stored energy fraction can no longer be unequivocally 
assigned to the field. 

Expression (3.2) is also interesting because it brings up the possibility of a 
negative total energy density being associated with a wave in a dispersive medium. 
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Although the idea of real electromagnetic fields being associated with a negative 
total energy density may seem unphysical, we will see in the next two sections that 
it has a clear and simple interpretation. 



3.2 Anomalous dispersion and temporary emis- 
sion 

Figure 3.1 depicts a planar pulse approaching and crossing an anomalously dispersive 
slab at normal incidence. For simplicity, we take the slab to have the properties 
a ~ 0, Cf ~ 1, fi r ~ 1, and n g « — 1 over the narrow bandwidth of the incident 
pulse. If we watch one particular segment of the pulse (for example that segment 
colored red in Figure 3.1) we see that when it crosses S\, the first dotted line, two 
copies of it appear downstream. One is beyond the dispersive medium and moving 
in the original direction of propagation. The other is in the anomalously dispersive 
medium and moves upstream. Eventually, the original upstream segment collides 
with its backward copy and the two dissappear, leaving only one downstream copy, 
which is now beyond the second dotted line, S , 2 . 

If we look at the total field energy associated with any given segment of the pulse, 
we see that it changes with time. Before the pulse segment reaches the first dotted 
line, it has an initial field energy, U$. As it crosses the dotted line, two copies are 
made. Counting the field energy associated with all three copies of the segment, we 
get a total field energy of 3U . As the most advanced copy crosses the second dotted 
line, the two lagging copies annihilate each other and the total field energy goes back 
to Uq. Conservation of energy requires that the extra 2U$ be accounted for by an 
equivalent energy deficit somewhere. Since the only other entity postulated besides 
the original pulse is the medium, it must come from the medium. The net effect is 
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that the medium donates energy to the field for a limited time and is later paid back. 
In Section 3.3, we will see that the medium is only paid back because the overall 
pulse has the narrow spectrum that we originally assigned to it. 

So far, we have found that the extra field energy that is present when a segment 
of the pulse is propagating through the anomalously dispersive medium suggests an 
energy deficit in the medium. We now seek to localize that energy deficit. In our 
model, energy is moved by electromagnetic fields, as represented by the Poynting 
vector. We can find the energy deficit associated with a particular slice of the medium 
by comparing the total energy that has entered it to the total energy that has left 
it. Because the nondispersive quantities associated with the medium have roughly 
their vacuum values, we know that the Poynting vector points downstream in the 
medium just as it does in the vacuum. 

Thus, to find the energy deficit associated with our slice, we compare the time 
integrated Poynting vector at the downstream surface of our slice, which corresponds 
to the energy which has left the slice, to the time integrated Poynting vector at the 
upstream surface of our slice, which corresponds to the energy which has entered the 
slice. The difference between these two quantities gives the total energy deficit in 
the medium between the two surfaces. 

The Poynting vector history associated with a surface in the medium is easy to 
visualize. If the surface were in the vacuum, the history would be given by the form 
of the pulse downstream from the surface. In our special medium, it corresponds to 
the form of the pulse that is upstream of the surface. In the geometry of Figure 3.1, 
the correspondence also ends at the line a, where annihilation destroys the record. 

When we subtract the energy that has crossed the first surface from the energy 
that has crossed the second surface, what we are left with corresponds to the energy 
between the two surfaces. In other words, the total energy deficit in our slice is equal 
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to the field energy in that slice. Applying this reasoning to a thin slice suggests that 
if the field energy density is Ufi e u and the density of energy absorbed by the medium 
is u exc hange, then we know for our scenario that 

"'field ' "•exchange "'field: 

Or tUat U exc y ian g e = ZU field- 

If rig = —2 (but all other quantities remain unaltered), the field energy relates 
to the Poynting vector history differently because a given portion of the pulse takes 
two times as long to cross a surface. Taking this into account, we find for this case 
that 

11 field i U eX change ^U field- 

Allowing n g to take an arbitrary value gives 

11 field i U eX change TlgU field- 

If we allow n to take values other than 1, we can apply this same logic. The longitu- 
dinal scaling of a pulse associated with n g must now be normalized by n and instead 
of u fi e id and u exc hange we can work in terms of a non-dispersive energy density (u ), 
given by Eq. 3.3, and a dispersive energy density (udi sp ), corresponding to the energy 
dispersively stored (or donated) by the medium. Then we can write 

n 
U + Udisp = —u - 
n 

Denoting the total energy associated with the propagating wave as u, we can rewrite 
this as 

n g 
u = —u , 
n 

which is exactly Eq. 3.2. 
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3.3 Temporary absorption and emission and the 
instantaneous spectrum in a lossless disper- 
sive medium 



When we calculate a pulse spectrum, we take an integral over the entire pulse. 
However, an atom responding to the pulse can only respond to that portion that has 
interacted with up to any particular point in time. Thus, the atom sees and interacts 
with an "instantaneous spectrum" that differs from the spectrum of the pulse as a 
whole. Peatross et al. put it this way [10]: 



The principle of causality requires a medium experiencing a pulse to be 
prepared for an abrupt termination of the field at any moment, in which 
case further exchange of energy with the field cannot take place. Such 
a termination produces a truncated waveform that generally contains a 
wider range of spectral components than are present in the pulse taken in 
its entirety. This momentary spectrum can lap onto nearby absorbing or 
amplifying resonances. The medium accordingly attenuates or amplifies 
this perceived spectrum. As the medium experiences the waveform, it 
continually reassesses the spectrum and thereby treats the front and the 
rear of the pulses differently. 



In this section we briefly review their results and seek to apply them to absorp- 
tionless dispersion in the narrow-band limit. 
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3.3.1 Instantaneous spectrum and energy exchange in a dis- 
persive medium 

Assuming, as we did in the previous section, that energy propagates only via the 
Poynting vector, we may write [10] 

V. B + £, (3.5) 

electromagnetic fields, where u, the total energy density, is given by 

u{t) = u fi e i d + u exchange + tt(-oo). (3.6) 

Here Uf ie \d represents the energy density of the electromagnetic fields, u exc h ange rep- 
resents the density of total energy absorbed or emitted by the medium, and u(— oo) 
represents the density of energy stored in the medium before the arrival of the field. 
Assuming that the medium is non-magnetic, Peatross et al. show, by imposing the 
requirement of causality on the linear electric susceptibility, that 



OG 



Uexchange = ^0 / \E t (u)\ UXi(f*>)du, (3.7) 

J —oo 

where 

1 /"* 

E t {u) = ^= dt'E{t')exp(iu)t'), (3.8) 

and E(u) and x(w) na d previously been defined in a way that matches the Fourier 
transform convention of Eq. (3.8). In other words, the u eX change is completely deter- 
mined by the interaction of the instantaneous power spectrum with the imaginary 
portion of the susceptibility. 

A narrow band pulse traveling through a dispersive, absorptionless medium has an 
instantaneous spectrum that narrows as the pulse amplitude decreases and broadens 
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as it increases. Although the medium is absorptionless with respect to the pulse 
as a whole, the fact that the medium is dispersive means in the spectral wings the 
medium is active (in the case of fast light), or absorbing (in the case of slow light). 

3.3.2 The instantaneous spectrum and narrow band pulses 
in dispersive absorptionless media 

In Section 3.1 we found a form for the total energy density associated with a pulse 
in a dispersive absorptionless medium and noted that the total energy density could 
be negative. In Section 3.2 we showed that the nature of the kinetic propagation of 
a pulse through an absorptionless anomalously dispersive medium suggests that the 
medium must temporarily donate energy to the electromagnetic field. We now see 
that the instantaneous spectrum provides the mechanism by which this occurs. 

An important point is that temporary absorption and temporary emission as dic- 
tated by the spontaneous spectrum are really just ordinary absorption and ordinary 
emission. The fact that they are temporary is specific to the way that later portions 
of the pulse interact with the medium and is guaranteed only when the overall pulse 
spectrum remains within the narrow bandwidth for which the absorption is close to 
zero. If a pulse is truncated, the energy that has been either absorbed or emitted 
remains either absorbed or emitted. 

This equivalence between dispersive energy exchange with the medium and ab- 
sorption and emission is useful in answering simple questions. For example, what 
kind of momentum exchange should we expect to see between field and medium 
when a pulse propagates through a dispersive medium? Because we know how mo- 
mentum exchange works under absorption and emission and because we know that 
light removed from the field or added to the field by dispersion has been temporarily 
absorbed or emitted, we can understand this momentum exchange in simple terms. 
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Chapter 4 discusses momentum in a dispersive medium in more detail. 
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n 9 =l 




Figure 3.1: A pulse interacts with an anomalously dispersive slab with a group 
index of —1 and with negligible absorption or gain. The two solid lines labeled a and 
c delimit changes in the group index. They also represent lines of symmetry. The 
magnitude of the Poynting vector is symmetric in form around line a between the 
dotted line Si and the solid line c. There is a similar symmetry about line c between 
lines a and 5*2. The portions of the pulse that are outside of the dotted lines have no 
copies. Those portions inside the dotted lines are represented in triplicate. Portions 
of the pulse between a and c move backwards. Watching the pulse propagate, it 
seems as if two pulses are created at the solid line labled c and then two pulses are 
later annihilated at the solid line labeled a. 
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Chapter 4 



Momentum in a dispersive medium 



4.1 Abstract 

When the effects of dispersion are included, neither the Abraham nor the Minkowski 
expression for electromagnetic momentum in a dielectric medium gives the correct 
recoil momentum for absorbers or emitters of radiation. The total momentum den- 
sity associated with a field in a dielectric medium has three contributions: (i) the 
Abraham momentum density of the field, (ii) the momentum density associated with 
the Abraham force, and (iii) a momentum density arising from the dispersive part 
of the response of the medium to the field, the latter having a form evidently first 
derived by D.F. Nelson [Phys. Rev. A44, 3985 (1991)]. All three contributions 
are required for momentum conservation in the recoil of an absorber or emitter in a 
dielectric medium. We consider the momentum exchanged and the force on a polar- 
izable particle (e.g., an atom or a small dielectric sphere) in a host dielectric when 
a pulse of light is incident upon it, including the dispersion of the dielectric medium 
as well as a dispersive component in the response of the particle to the field. The 
force can be greatly increased in slow-light dielectric media. 
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4.2 Introduction 

Electromagnetic momentum in a dielectric medium is a subject with a very extensive 
literature, especially in connection with its different formulations. The two most 
favored forms by far are those of Abraham and Minkowski; as aptly remarked in 
a recent paper [35], "There is ... a bewildering array of experimental studies and 
associated theoretical analyses which appear to favor one or other of these momenta 
or, indeed, others." An aspect of this subject that has received surprisingly little 
attention concerns the effects of dispersion on the Minkowski and Abraham momenta 
and on the electromagnetic forces on polarizable particles. The intent of the present 
paper is to address such effects, which might help to clarify the physical interpretation 
of the Abraham and Minkowski momenta and the distinction between them. 

We first review briefly the Abraham and Minkowski momenta for the situa- 
tion usually considered — a dielectric medium assumed to be dispersionless and non- 
absorbing at a frequency u. The Abraham and Minkowski momentum densities are 
respectively 

P A = — E x H and P M = DxB (4.1) 

c 2 

in the standard notation for the fields on the right-hand sides. We will take the 

permeability u. to be equal to its vacuum value /io, which is generally an excellent 

approximation at optical frequencies. For single photons the magnitudes of the 

Abraham and Minkowski momenta are given by (see Section II) 

1 hu fiuj 

Pa = and p M = n — , (4.2) 

n c c 

where n is the refractive index at frequency u. From D = eo^ 2 E it follows that 

dP M dP A+fA ^ (43) 



(4.4) 



dt 




dt 


1 ■"■ ! 




where 










f A = 


1 

7- 


(n 2 - 
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is the Abraham force density. For single-photon fields the momentum p A associated 
with the Abraham force is [{n 2 — l)/n]fouj/c, and (4.3) becomes pu = Pa + P A - 

The Abraham momentum is generally regarded as the correct momentum of the 
electromagnetic field [3G] , whereas the Minkowski momentum evidently includes the 
momentum of the dielectric medium as well as that of the field. Ginzburg [37] calls 
Pm the momentum of a "photon in a medium," and notes that its use, together 
with energy and momentum conservation laws, yields correct results for Cerenkov 
radiation as well as the Doppler shift. Experiments appear by and large to indicate 
that it is the momentum nhu/c per photon that provides the recoil and radiation 
pressure experienced by an object immersed in a dielectric medium [3 ]. However, 
when dispersion (dn/du) is accounted for, nhu/c is not the Minkowski momentum 
of a photon, as we review in the following section. 

This paper is organized as follows. In the following section we briefly discuss 
the generalization of the Abraham and Minkowski momenta to the case of a disper- 
sive dielectric medium [. ] and consider two examples: (i) the Doppler shift in a 
dielectric medium [40] and (ii) the displacement of a dielectric block on a friction- 
less surface due to the passage of a single-photon field through it [ ] . A consistent 
description of momentum transfer in these examples requires that we account for 
momentum imparted to the medium. In Section 4.4 we calculate the force exerted 
by a quasimonochromatic plane wave on a polarizable particle and on a dispersive 
dielectric medium modeled as a continuum, and obtain a dispersive contribution to 
the latter in agreement with an expression that, to the best of our knowledge, was 
first derived, in a rather different way, by Nelson [42]. In Section 4.5 we consider the 
momentum exchange between a plane-wave pulse and an electrically polarizable par- 
ticle immersed in a nonabsorbing dielectric medium, and show that this momentum 
depends on both the dispersion of the medium and the variation with frequency of 
the polarizability; in particular, in slow-light media it can be large and in the direc- 
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tion opposite to that in which the field propagates. Section 4.6 presents derivations 
of some results relevant to Section 4.7, where we generalize the results of Section 4.5 
to include absorption and discuss the forces exerted by a pulse on a small dielectric 
sphere in a host slow-light medium. Section 4.8 briefly summarizes our conclusions. 

4.3 Abraham and Minkowski Momenta for Dis- 
persive Media 

We first recall the expression for the total cycle-averaged energy density when a plane- 
wave monochromatic field [E = E^e - '"*, H = H w e~ iaji , H 2 = (e//i )E 2 ] propagates 
in a dispersive dielectric at a frequency u at which absorption is negligible [43]: 
1 



U = l 



^-(eu)El + f , Hl 



(4.5) 



or equivalently, in terms of E^ and the group index n g = d(noj)/duj, 

u = -e nn g El (4.6) 

When the field is quantized in a volume V, u is in effect replaced by qhw/V, where q 

is the expectation value of the photon number in the volume V\ therefore, from (4.6), 

E^ is effectively 2hw / (eonn g V) per photon. Thus, for single photons, the Abraham 

momentum defined by (4.1) is 

n 1 2tvuj 1 tko . , „. 

^ = -c2 t "7^^ V -n g -- (47) 

Similarly, 

n 2 Hcu , , . 

Pm = , 4.8 

n g c 

which follows from the definition in (4.1) and the relation D = eo«- 2 E; thus Pm = 
n 2 pA- These same expressions for pa and pu can of course be obtained more formally 
by quantizing the fields E, D, H, and B in a dispersive medium [3 
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Two examples serve to clarify the differences among the momenta involved in the 
momentum exchange between light and matter. The first example is based on an 
argument of Fermi's that the Doppler effect is a consequence of this momentum ex- 
change [40], as follows. Consider an atom of mass M inside a host dielectric medium 
with refractive index n(u). The atom has a sharply defined transition frequency 
uq and is initially moving with velocity v away from a source of light of frequency 
u. Because the light in the atom's reference frame has a Doppler-shifted frequency 
u(l — nv/c) determined by the phase velocity (c/n) of light in the medium, the atom 
can absorb a photon if u(l — nv/c) = uq, or if 

u^u (l + nv/c). (4.9) 

We denote the momentum associated with a photon in the medium by p and consider 
the implications of (nonrelativistic) energy and momentum conservation. The initial 
energy is Ei = hu + |Mt> 2 , and the final energy, after the atom has absorbed a 
photon, is |Mi/ 2 + Hcjq, where v' is the velocity of the atom after absorption. The 
initial momentum is p + Mv, and the final momentum is just Mv '. Therefore 

X -M{v' 2 - v 2 ) S Mv(v' -v) = Mv(p/M) = h(cu - w ), (4.10) 

or u = w + pv/H. From (4.9) and u = uo we conclude that 

p = n^. (4.11) 

Thus, once we accept the fact that the Doppler shift depends on the refractive 
index of the medium according to Eq. (4.9), we are led by energy and momentum 
conservation to conclude that an atom in the medium must recoil with momentum 
(4.11) when it absorbs (or emits) a photon of energy fi/jj. Momentum conservation 
in this example is discussed in more detail below. 

In our second example we consider, following Balazs [41], a rigid block of mass M, 
refractive index n, and length a, initially sitting at rest on a frictionless surface. A 
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single-photon pulse of frequency u passes through the block, which is assumed to be 
nonabsorbing at frequency u and to have anti-reflection coatings on its front and back 
surfaces. The length a of the block is presumed to be much larger than the length of 
the pulse. If the photon momentum is p in inside the block and p out outside, the block 
picks up a momentum MV = p out — p\ n when the pulse enters. If the space outside 
the block is vacuum, p out = mc, where m = E/c 2 = ftw/c 2 . Similarly p- in = mv p , 
where v p is the velocity of light in the block. Without dispersion, v p = c/n and the 
momentum of the photon in the block is evidently p- m = mc/n = hu/nc. The effect 
of dispersion is to replace v p = c/n by v g = c/n g and p- m = hw/nc by p in = huj/n g c. 
With or without dispersion, this example suggests that the photon momentum in 
the medium has the Abraham form. Note that the essential feature of Balazs's 
argument is simply that the velocity of light in the medium is v p (or, more generally, 
v g ). This, together with momentum conservation, is what leads him to conclude that 
the momentum of the field has the Abraham form. 

This prediction can in principle be tested experimentally. Conservation of mo- 
mentum requires, according to Balazs's argument, that MV = m(c — v g ). When 
the pulse exits the block, the block recoils and comes to rest, and is left with a net 
displacement 



m a hw 

M {C ~ Vg) v~ = M? 



Ax = VAt = -(c - v g )- = —{n g - l)a (4.12) 



as a result of the light having passed through it. This is the prediction for the net 
displacement based on the momentum pa given in (4.7). If the photon momentum 
inside the block were assumed to have the Minkowski form n 2 hu/cn g given in (4.8), 
however, the displacement of the block would in similar fashion be predicted to be 

Ax=— a(n 9 -n 2 ), (4.13) 

and if it were assumed to be nhu/c, as in Eq. (4.11), the prediction would be that 
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the net displacement of the block is 
Mc 2 



Ax = -j^-^arig^l - n). (4.14) 



These different assumptions about the photon momentum can lead to different pre- 
dictions not only for the magnitude of the block displacement but also for its direc- 
tion. 

The first (Doppler) example suggests at first thought that the momentum of 
the photon is nhid/c [Eq. (4.11)], while the second (Balazs) example indicates that 
it is hw/rigC. Let us consider more carefully the first example. There is ample 
experimental evidence that the Doppler shift is nvu/c regardless of dispersion, as we 
have assumed, but does this imply that the momentum of a photon in a dielectric 
is in fact nhw/c! We will show in the following section that the forces exerted by 
a plane monochromatic wave on the polarizable particles of a dielectric result in a 
momentum density of magnitude 

Pmed = ^n(nn g - 1)E* = (n - ~) — y\ (4-15) 

the second equality applies to a single photon, and follows from the replacement 
of E^, by 2hw/(e nn g V), as discussed earlier. Now from the fact that the Doppler 
shift implies that an absorber (or emitter) inside a dielectric recoils with momentum 
nhw/c, all we can safely conclude from momentum conservation is that a momentum 
nhw/c is taken from (or given to) the combined system of field and dielectric. Given 
that the medium has a momentum density (4.15) due to the force exerted on it by 
the propagating field, we can attribute to the field (by conservation of momentum) 
a momentum density 

fhd 1 1 hid 1 
n — --Pmed = 77= PA- (4.16) 

c V n g c V 

That is, the momentum of the field in this interpretation is given by the Abraham 
formula, consistent with the conclusion of the Balazs thought experiment. The re- 
coil momentum nhu/c, which in general differs from both the Abraham and the 
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Minkowski momenta, evidently gives the momentum not of the field as such but of 
the combined system of field plus dielectric. It is the momentum density equal to 
the total energy density u = Hoo/V for a monochromatic field divided by the phase 
velocity c/n of the propagating wave. As already mentioned, experiments on the 
recoil of objects immersed in dielectric media have generally indicated that the recoil 
momentum is nhu/c per unit of energy Hu of the field, just as in the Doppler effect. 
But this should not be taken to mean that nhu/c is the momentum of a "photon" 
existing independently of the medium in which the field propagates. Regardless of 
how this momentum is apportioned between the field and the medium in which it 
propagates, the important thing for the theory, of course, is that it correctly predicts 
the observable forces exerted by electromagnetic fields. We next turn our attention 
specifically to the forces acting on polarizable particles in applied electromagnetic 
fields. 



4.4 Momenta and Forces on Polarizable Particles 

We will make the electric dipole approximation and consider field frequencies such 
that absorption is negligible. Then the induced electric dipole moment of a particle 
in a field of frequency u is d = a(w)E aj exp(— iut), and the polarizability a(u) may 
be taken to be real. With these assumptions we now consider the forces acting on 
such particles in applied, quasi-monochromatic fields. 

We begin with the Lorentz force on an electric dipole moment d in an electro- 
magnetic field [44]: 

F = (d-V)E + dxB 

= (d-V)E + dx(VxE) + ^(dxB) 

at 

= F B + F B , (4.17) 
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where we define 



F E = (d-V)E + dx (VxE) 



(4.18) 



F fl = £(dxB). 



(4.19) 



In writing the second equality in (4.17) we have used the Maxwell equation dB/dt = 
— V x E. The dipole moment of interest here is induced by the electric field. Writing 



E = £o(r,t)e 



-iujt —iwt 



dA£ (r,A)e 



-iAt 



(4.20) 



in which \dSo/dt\ <C u)\£q\ for a quasi-monochromatic field, we approximate d as 
follows: 



d(r,£) 



dAa(u + A)£ {r,A)e~ i{ul+A)t 
dA[a(u) + Aa'(u)]S (r, A) e - i(w+A)i 



a(u)£o(r, t) + ia'(u 



d£o 
dt 



-iujt 



(4.21) 



Here a' = da/du and we assume that higher-order dispersion is sufficiently weak 
that terms d m a/du m can be neglected for m > 2. Putting (4.21) into (4.18), we 
obtain after some straightforward manipulations and cycle-averaging the force 
'1 



-a(U))\S\ 



4 



1 1 Ml Cl2 



(4.22) 



where £ and k are defined by writing £ (r,t) = £(r,t)e l]i ' v . Since the refractive 
index n of a medium in which local field corrections are negligible is given in terms 
of a by n 2 — 1 = Na/eo, N being the density of dipoles in the dielectric, we have 
a' = (2n€o/N)(dn/doj) and 
'1 



-a(u})\£\ 



+ iL kn ^^.|P|2 

+ 2N KU dudt ltl • 



(4.23) 
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The first term is the "dipole force" associated with the energy W = — \a{uj)E 2 
involved in inducing an electric dipole moment in an electric field: 

W = - d-dE = -a(co) E • dE = — a(u)E 2 . (4.24) 

Jo Jo 2 

The second term in (4.23) is nonvanishing only because of dispersion (dn/du ^ 0). 

It is in the direction of propagation of the field, and implies for a uniform density N 

of atoms per unit volume a momentum density of magnitude 

Po = \ton 2 ~\£\ 2 = ^n\n g - n)\S\ 2 , (4.25) 

Z CLUJ C Z C 

since k = n(u)u/c. This momentum density comes specifically from the dispersion 
(dn/du) of the medium. 

The force F# defined by (4.19), similarly, implies a momentum density P A im- 
parted to the medium: 

P^iVdx B. (4.26) 

As the notation suggests, this momentum density is associated with the Abraham 
force density (4.4). The result of a straightforward evaluation of P A based on (4.21) 
and V x E = -dB/dt is 

P A = ie (n 2 - 1)-|£| 2 , P A = \ C Mn 2 ~ 1)|£| 2 , (4.27) 

Z Ul z c 

when we use k • E = and our assumption that |£o| "C w|£o|- The magnitude of the 
total momentum density in the medium due to the force of the field on the dipoles 
is therefore 



3 med = Pd + P A = f c [n 2 {n g -n)+ n{n 2 - 1)] \£\ 



£2. r^2 
2c 

= e fn{nn g - l)\£\ 2 (4.28) 

in the approximation in which the field is sufficiently uniform that we can ignore the 
dipole force V[t«|£| 2 ]. 
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The complete momentum density for the field and the medium is obtained by 
adding to (4.28) the Abraham momentum density Pa of the field. According to 
(4.1), Pa = (eo/2c)n|£| 2 , and so the total momentum density is 

P A + P D + P A = € f c [n + n(nn g - 1)}\£\ 2 = ^n 2 n g \E\ 2 (4.29) 

if the dipole force is negligible. To express these results in terms of single photons, 
we again replace |£ | 2 by 2hu/(e nn g V); then (4.29) takes the form 

PA + PD+p A = n—- (4.30) 

c V 

consistent with the discussion in the preceding section. This is the total momentum 
density per photon, assuming that the dipole force is negligible. The momentum 
density of the medium per photon follows from (4.28): 

, a e o / , 2hw 

Pmed = Pd+P =—n{nn g -l) — 

2c y nn g e V 

, l , hw l . , „ . 

= (n--) — - } 4.31 

rig c V 

as stated earlier [Eq. (4.15)]. 

Consider the example of spontaneous emission by a guest atom in a host dielectric 
medium. The atom loses energy fkvo, and the quantum (photon in the medium) 
of excitation carries away from the atom not only this energy but also a linear 
momentum nhw/c [Eq. (4.30)]. The atom therefore recoils with momentum nfiw/c 
[45]. 

The momentum density (4.25) was obtained by Nelson [ 12] in a rigorous treat- 
ment of a deformable dielectric based on a Lagrangian formulation; in the present 
paper a dielectric medium is treated as an idealized rigid body. From a microscopic 
perspective, this part of the momentum density of the medium is attributable di- 
rectly to the second term on the right-hand side of (4.21), i.e., to the part of the 
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induced dipole moment that arises from dispersion. In the Appendix the relation of 
this term to the formula (4.5) for the total energy density is reviewed; the term is 
obviously a general property of induced dipole moments in applied fields. Consider, 
for example, a two-level atom driven by a quasi-monochromatic field with frequency 
u far-detuned from the atom's resonance frequency uq. In the standard u, v no- 
tation for the off-diagonal components of the density matrix in the rotating-wave 
approximation [46], 

u(t)- i v(t) = ± X (t) + ^ + ..., (4.32) 

where x{t) is the Rabi frequency and A is the detuning. The polarizability is pro- 
portional to 1/A in this approximation, and therefore (4.32) is just a special case of 
(4.21). 



4.5 Momentum Exchange between a Light Pulse 
and an Induced Dipole 

We next consider the momentum exchange between a plane-wave pulse and a single 
polarizable particle. We will assume again that the particle is characterized by a real 
polarizability a(u) and that it is surrounded by a host medium with refractive index 
rib(oj). The electric field is assumed to be 

E(z, t) = S(t - z/v bg ) cos(ut - kz), (4.33) 

with k = rib(u))u)/c and group velocity Vb g = c/nbg, rib g = (d/du)(ujnb). 



The force acting on the particle is F^ + F^. F# reduces to ^a(u)(k/u)(d/dt)\£\ 



2 



obtained by multiplying (4.27) by a volume V describing the pulse, replacing n 2 — 1 
by Na/eo with NV = 1 for the single particle, and differentiation with respect to 



49 



Chapter 4. Momentum in a dispersive medium 

time. F# follows from (4.22). Then the force acting on the particle is in the z 
direction and has the (cycle-averaged) magnitude 

+ LaMn„(u>)^e 2 , (4.34) 

where now we retain the dipole force, given by the first term on the right-hand side. 
The momentum of the particle at z at time T is 

T 



f 1 f d 

p = Fdt = -a / —£ 2 (t — z/vb g )dt 

J — oo ■* J— oo ^-^ 



1 /" T (9 

+— -a'n&o; / ^ 2 ( t ~ z l v b g )dt 

i /-^a 



+ Y aHh \ ~ft£ (* _ z l v bg)dt 



■X) 



— a — £ H awe + -a — £ 

4 v bg 4c 2 c 

— [(2n 6 - n bg )a + n^a^^T - z/v bg ). (4.35) 



Hinds and Barnett [ ] have considered the force on a two-level atom due to a 
pulse of light in free space. In this case n b = n bg = 1 and (4.35) reduces to 

p = —[a + ua'}£ 2 . (4.36) 

Following Hinds and Barnett, we argue that a pulse occupying the volume V in the 
neighborhood of the atom in free space corresponds to a number q = ^e^V/hw of 
photons, so that 

p=— [a+wa']— g. (4.37) 

a = eo(n 2 — l)/N, where n is the refractive index in the case of N polarizable particles 
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per unit volume. Then 

£o(n 2 — 1) 2eon dn 



1 
V z 



H —u- 

N N du 



frw 

— q 

c 



2c ^ 

r dn.hu: Tr fuxi , 

n - 1 + uo—} — q = K — q. 4.38 

du c c 



This is the momentum imparted to the particle, which implies a change in field 
momentum per photon equal to 

fuxi r rr , hw l tvjj , 

— [l-K} = —- - = 4.39 

c L J c 1 + K n g c v ' 

\i\K\ <^i 1, where n g = (d/du)(nu). As in the case of a two-level atom considered by 
Hinds and Barnett, this corresponds to the Abraham momentum; our result simply 
generalizes theirs in replacing n by n g in the expression for the change in photon 
momentum. 

In the case of a polarizable particle in a host dielectric rather than in free space 
we obtain, from (4.35), 

p=JL[(2-^) a + U0 {\, (4.40) 

2e c 2 n b 

where the intensity / = (l/2)ceoni,£ 2 . If dispersion in the medium and in the po- 
larizability of the guest particle are negligible, we can set n^g = n and a' = 0, and 
then (4.40) reduces to a well known expression [47]. However, this momentum can 
be large in a slow-light medium (n& 9 large), for example, because the gradient of 
the field (4.33) responsible for the dipole force on the particle is large [48]; this is a 
consequence of the spatial compression of a pulse in a slow-light medium. We discuss 
this case further in Section 4.7. But first we return to some other well known results 
that are relevant there. 
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4.6 Electric Dipole Radiation Rate and Rayleigh 
Scattering 

A Hertz vector II (r, u) can be defined for a dielectric medium, analogous to the case 
of free space [49] , by writing the electric and magnetic field components at frequency 

uj as 

E(r, w) = A; 2 [e 6 (u;)/eo]II(r, u) + V[V ■ II(r, w)], (4.41) 

H(r, u) = -iue b {u)V x n(r, u). (4.42) 

Here k = u/c and we denote by e b (u) the (real) permittivity of the dielectric. We 
will be interested here in a dipole source inside the "background" dielectric medium. 
The identifications (4.41) and (4.42) are consistent with the propagation of a wave 



of frequency u with the phase velocity c/n b (u) in the medium [n b (u) = \/e&(w)/eo], 
as will be clear in the following. 

The curl of E(r,w) in (4.41) is simply 

V x E(r, u) = k 2 [e b (uj)/eo}V x n(r, w), (4.43) 

since the curl of a gradient is zero. Now apply the curl operation to this equation, 
assuming no free currents and therefore V x H(r,w) = -iwD(r,w): 

V x (V x E) = iLOj2 V x H = uj 2 jj, D 

= k 2 [e b (co)/e ]V x (V x n) 

= fc 2 [e b (o;)/eo][V(V-n)-V 2 n], (4.44) 

implying 

V 2 n = -^/ioD + V(V ■ n) = --D + [E - -fc 2 n], (4.45) 

e& k e b e 
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V 2 n + k 2 II = E - -D, k 2 = kle b {u)/e = n 2 h (u)u 2 /c 2 . (4.46) 

If D(r,w) = e&(u;)E(r, u), the right-hand side is zero, and all we have done is red- 
erived what we already know: the field propagates with phase velocity u/k(u) = 
c/n^u). Suppose, however, that within the medium there is a localized source char- 
acterized by a dipole moment density P s (r,w) = po(u)5 3 (r). Then D = e^E + P s 
and 

V 2 n + k 2 n = -- P oH5 3 (r). (4.47) 

The solution of this equation for II(r, u) is simply 

and from this one obtains the electric and magnetic fields due to the source in the 
medium. In the far field, assuming po = pz and letting 9 be the angle between the 
z axis and the observation point, 

E e = p?-sme — , (4.49) 

47re r 



n b k 2 p [1^ . e ikr 

Ha > = — — — ,/ — sinS , (4.50) 

47re V A*o r 

in spherical coordinates. The Poynting vector S = E x H implies the radiation rate 

P = ^4. (4.5D 

lzTreocr 

analogous to the fact that the spontaneous emission rate of an atom in a dielectric 
without local field corrections is proportional to the (real) refractive index at the 
emission frequency. 
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Polarizability of a Dielectric Sphere 

Suppose, somewhat more generally, that the source within the medium occupies a vol- 
ume V and is characterized by a permittivity e s (u). Then D(r, oS) = e(r, w)E(r, w), 
where e = e s (u) within the volume V occupied by the source and e(r, u) = e b (u) 
outside this volume, and 



V 2 n + k 2 U = [1 - e(r, w)/c 6 (w)]E. 



(4.52) 



The solution of this equation is 



U(r,u) 



1 

All 



1 



e. w 



e&(w) 



ik\r-r'\ 

d 3 r'E(r',u)r r 

v ; |r-r'| 



(4.53) 



Suppose further that the extent of the volume V is sufficiently small compared to a 
wavelength that we can approximate (4.53) by 



n(r,w) = - 



47T 



1- 



e&(w) 



Akr 



VE iQS (u)- 



(4.54) 



with r the distance from the center of the source (at r = 0) to the observation point 
and Ei ns (w) the (approximately constant) electric field in the source volume V. This 
has the same form as (4.48) with Po(oj) = €b(uj)[e s (u) / e b {u) — l]VE ins (w). In other 
words, II(r, u) has the same form as the Hertz vector for an electric dipole moment 



p (w) = [e,(w) - e b (u)]VE- ms (uj). 



(4.55) 



Consider, for example, a small dielectric sphere of radius a: V = 4ira 3 /3. The field 
inside such a sphere is E ins (u;) = [3e b /(e s + 2e b )]E b (uj) , where E b (oj) is the (uniform) 
electric field in the medium in the absence of the source. The dipole moment (4.55) 
in this case is therefore related to the external field E out (w) by Po(u>) = a(w)E out (w), 
where the polarizability 



a(u) = 4:iie b 



e g -e b 
e s + 2e fe 



(4.56) 
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Rayleigh Attenuation Coefficient 

The cross section for Rayleigh scattering for an ideal gas of refractive index n{uS) 
can be deduced as follows [ ] . An electric field E cos ut induces an electric dipole 
moment p(t) = a(u)E cos cut in each of N isotropic, polarizable particles per unit 
volume, each particle having a spatial extent small compared to a wavelength. The 
power radiated by this dipole is, from Eq. (4.51), 

(IU ''"" n(u)-^-y(uj)El = a R (co)I, (4.57) 



dt 127re c 3 

where W ra d denotes energy of the radiated field, I = ^n(oj)ce El is the intensity of 
the field incident on the dipole, and 

c,rM = ^(^)VM-1] 2 (4-58) 

is the (Rayleigh) scattering cross section. We have assumed that local field correc- 
tions are negligible and used the formula n 2 (oj) — 1 = iVa(w)/e to express a R (u) in 
terms of the refractive index n(u>). The attenuation coefficient is then 

a R = Na R =^Q\n^)-lf. (4.59) 

Rosenfeld [51] obtains instead 

because he does not account for the factor n{uS) in the dipole radiation rate (4.57). 
Rayleigh's derivation of (4.59) follows essentially the one just given, but the factor 
n(uj) appears in neither the dipole radiation rate nor the expression for the intensity 
(or actually, in his derivation, the energy density) [ .]. In practice the difference 
between (4.59) and (4.60) is negligible for the case assumed here of a dilute medium 
[53]. 
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4.7 Force on a Dielectric Sphere 

The expression (4.34) for the force on a polarizable particle in a field (4.33) may be 
generalized to allow for absorption by the particle simply by taking the polarizability 
a(u>) in (4.21) to be complex. Assuming again that 8 is slowly varying in time 
compared to exp(-iut) , and slowly varying in space compared to exp(ikz), we obtain 

F = — [(2n b - n bg )a R + n b ua' R ] — \£\ 2 + -n b -a/|£| 2 , (4.61) 

where r = t — n bg z/c and an and olj are the real and imaginary parts, respectively, of 
a{yS). If we replace n bg by n b and take a' R = 0, we recover results that may be found 
in many previous works when absorption is assumed to be negligible [47]. The last 
term in (4.61) is the absorptive contribution to equation (7) of a paper by Chaumet 
and Nieto-Vesperinas [54] when the field is assumed to have the form (4.33). 

The polarizability in the case of a dielectric sphere of radius a much smaller than 
the wavelength of the field is given by (4.56). Dispersion affects the force (4.61) both 
through the group index (n bg ) of the host dielectric medium and the variation of the 
real part of the sphere's polarizability with frequency (a' R ). The latter depends on 
both the intrinsic frequency dependence of the permittivity of the material of the 
sphere and the frequency dependence of the refractive index of the host medium. If 
these dispersive contributions to the force exceed the remaining two contributions to 
the force (4.56), 

F - h {~ aRrib ° + nbUja '^ ^ |£|2 - (4 - 62) 



Using (4.56) for this case, we obtain 

-» ^..;''"L 2 ^| 2 (4.63) 



37reoa 3 n2 s n b 9 

<- " (n 2 s + 2nl) 2 dr ' 



if the dispersion of the dielectric material constituting the sphere is much smaller 
than that of the host dielectric medium, i.e., if de s /du> ^C de b /du. (Here n s is the 
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refractive index at frequency u of the material of the sphere.) This result implies 
that, in the case of a slow-light host medium {n bg 3> 1), the force on the sphere can 
be much larger than would be the case in a "normally dispersive" medium, and is in 
the direction opposite to that in which the field propagates. 

The simple formula (4.63), and similar expressions obtained in other limiting 

cases of (4.61), obviously allow for a wide range of forces when a pulse of radiation 

is incident on a dielectric sphere in a host dielectric medium. Here we make only a 

few remarks concerning the last term in (4.61). Although we have associated this 

contribution to the force with absorption, such a force appears even if the sphere 

does not absorb any radiation of frequency u. This is because there must be an 

imaginary part of the polarizability simply because the sphere scatters radiation and 

thereby takes energy out of the incident field. According to the optical theorem in 

this case of scattering by a nonabsorbing polarizable particle that is small compared 

to the wavelength of the field, the imaginary part of the polarizability is related to 

the complete (complex) polarizability as follows [55]: 

1 2c«j 3 
a/H = —^n b \a(u)\ 2 . (4.64) 

Then the force proportional to a,i(oS) in (4.61) is 

^ 1 kW . _ 9 87r /u\ 4 n^I ( e s — e b \ fi . 

F scat = -n 5 b - ai \£\ 2 = — (-) -M s — M a 6 , (4.65) 



2 c 3 \cJ c \e s + 2eft / 

which is just the well known "scattering force" [56] on a dielectric sphere in a medium 
with refractive index n&, which may be taken to be real in the approximation in which 
the field is far from any absorption resonances of the sphere. 



4.8 Conclusions 

In this attempt to better understand the different electromagnetic momenta and the 
forces on electrically polarizable particles in dispersive dielectric media, we have made 
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several simplifications, including the neglect of any surface effects, the treatment of 
the medium as a nondeformable body, and the approximation of plane-wave fields. 
We have shown that conservation of momentum, even in seemingly simple examples 
such as the Doppler effect, generally requires consideration not only of the Abraham 
momentum and the Abraham force, but also of a contribution to the momentum 
of the medium due specifically to the dispersive nature of the medium. We have 
generalized some well known expressions for the forces on particles immersed in 
a dielectric medium to include dispersion. While we have presented arguments in 
favor of the interpretation of the Abraham momentum as the momentum of the field, 
our simplified analyses lead us to the conclusion that neither the Abraham nor the 
Minkowski expressions for momentum give the recoil momentum of a particle in a 
dispersive dielectric medium. Finally we have shown that the force exerted on a 
particle in a strongly dispersive medium is approximately proportional to the group 
index n& s , and can therefore become very large in a slow- light medium. 
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Nonstationary electromagnetics of 
controllably dispersive media 



Abstract 



Recent experiments have demonstrated that it is possible to alter the dispersion of a 
medium without significantly altering its absorption or refractive index and that this 
may be done while a wave propagates through the medium. This possibility opens up 
a new set of potential experiments to the field of nonstationary optics. We consider 
the basic kinetics of waves propagating through a medium whose group and phase 
velocities are a function of position and time. We compare the dynamics of waves 
propagating through two homogeneous media, one a nondispersive medium with a 
time dependent phase velocity, and one a dispersive medium with a time dependent 
group velocity and show that new dynamic effects accompany new kinetic ones. 
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5.1 Introduction 

The terms nonstationary, inhomogeneous, and dispersive describe media whose prop- 
erties vary with time, position, and frequency, respectively. In this paper, we consider 
transformations wrought upon waves as they propagate through nonstationary, inho- 
mogeneous media with variable dispersion. Our interest in this topic is motivated by 
experimental advances in the manipulation of dispersion through nonlinear optics. 
Recent experiments have demonstrated that it is possible to achieve either strong nor- 
mal dispersion (in the case "slow-light media [7, 9]") or strong anomalous dispersion 
(in the case of "superluminal media [ ]") in nearly transparent media using nonlin- 
ear optical effects. In either case, the strength of the dispersion, as experienced by a 
probe beam, depends on the intensity of one or more controllable auxiliary beams. 
If auxiliary intensities are changed while the probe beam is in transit, then the dis- 
persion, as perceived by the probe beam, is time dependent. These media, which we 
will refer to collectively as controllably dispersive media, thus form a newly accessible 
class of nonstationary media. Although the nonstationarity of a slow light medium 
has been used experimentally to adiabatically transform a "slow light" pulse to a 
"stopped" one [9, 57], an explicit connection between this application and the field 
of nonstationary electromagnetics has not, to our knowledge, been made. 

The field of nonstationary electromagnetics, although much less developed than 
its inhomogeneous counterpart, is now a venerable one. In 1958, Morgenthaler intro- 
duced propagation equations for electromagnetic waves in an isotropic nondispersive 
medium whose permittivity and permeability were allowed to vary with time (but 
not with space) [12]. Although the propagation equations are not generally analyti- 
cally soluble, he found useful solutions assuming, on the one hand, a step functional 
time dependence and, on the other hand, an adiabatic time dependence. He found, 
among other things, that the frequency of a wave would vary with the permittivity 
and permeability in such a way that wavelength would be preserved. 
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Over the next two decades, much of the progress of the field was in the Soviet 
Union, with a particular focus on the behavior of light in plasmas. We refer the 
reader to Stepanov and colleagues [58, 59] for an introduction to this literature and 
highlight here two aspects of it. First, the two simplified cases (step function and 
adiabatic changes) used by Morgenthaler continued to be used extensively, and were 
generalized from a purely nonstationary picture to one that allowed inhomogeneities 
to propagate at a fixed velocity according to the traveling wave law (/ = f(x — vt)) 
[58, 60]. This picture allows a single framework to unite the purely nonstationary 
case (v — >■ oo) 1 , a 1-dimensional purely inhomogeneous case (v — 0), and the case 
of an inhomogeneity (for example, an ionization front) moving at any velocity in 
between these two extremes. Second, a distinction between "kinetic" and "dynamic" 
phenomena was found to be useful in differentiating between those results of the 
theory that are general to all linear wave phenomena and those that are specific to 
particular media [60, 58, 61]. We employ both of these concepts in this paper. 

The advent of the laser, and particularly the pulsed laser, introduced a new way 
to dynamically alter the characteristics of a medium. Although they were apparently 
unaware of previous work in nonstationary electromagnetics, a few authors noted the 
possibility of using light to alter the phase velocity of propagating electromagnetic 
waves. In 1977, Lampe, Ott, and Walker noted that an ionizing laser might be swept 
across a gas to create a superluminal ionization front that could interact with a mi- 
crowave pulse in a nonstationary fashion [62]. In 1988, Wilks, Dawson, and Mori 
examined the problem of a wave propagating through a medium which is then quickly 
ionized by a high-intensity ultrashort pulse [ :]. Large frequency shifts were soon 
realized for waves reflecting off of a relativistically propagating ionization front [64] . 
These works founded a sub-literature surrounding nonstationary effects in plasmas 



1 Here we note a basic characteristic of nonstationary media: there is no maximum speed 
because changes in a medium at two different positions need not be causally connected to 
each other. 
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[65, 66, 67, 68, 69, 11]. We point the reader to Shvartsburg, who reviewed some of 
the recent work on nonstationary effects in plasmas [69]. We also highlight the work 
of Avitzour and Shvets, whose 2008 paper relies on controlled dispersion in nonsta- 
tionary media. They proposed a method for altering the spectral width of a pulse 
without changing its central frequency using controlled dispersion in a nonstationary 
magnetized plasma [ ]. We will look at this method from a kinetic perspective in 
the next section. 

In parallel with the more recent work on nonstationary plasmas, several authors 
have explored electromagnetic propagation through nonstationary media abstractly, 
in the tradition of Morgenthaler [70, 71, 7 ]. Dodonov, Klimov, and Nikonov studied 
the quantization of a linear, nonstationary medium and were able to quantitatively 
relate temporal changes in dielectric permittivity to photon generation [' ]. Bian- 
calana, Amann, Uskov, and Oreilly treated a 1-dimensional nondispersive nonsta- 
tionary medium by introducing a transmission matrix for moving interfaces. Using 
this matrix they generalized Bragg reflection to propagating interfaces and showed 
that temporal periodicity leads to k-vector bandgaps just as spatial periodicity leads 
to frequency bandgaps [ ]. Budko has found an analogy between a nonstationary 
nondispersive medium and an expanding Universe [72]. 

One surprising fact about recent works, both those focusing on nonstationary 
plasmas [62, 63, 64, 65, 66, 67, 68, 69, 11] and those dealing with abstract, nondisper- 
sive, linear media [70, 71, 72], is that they make no reference to the substantial Soviet 
literature from the 1960s and 1970s. One consequence of the disjointedness of which 
this fact is symptomatic is that basic principles have had to be rediscovered, multiple 
times, in incidental and typically incomplete ways. For example, a basic principle 
of wave propagation in nonstationary linear media is the importance and general- 
ity of kinetic effects. Simple spatial and temporal symmetries play an important 
role in determining the effects of nonstationary propagating waves. Aspects of these 
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symmetries have been rediscovered and used in many contexts [12, 60, 63, 71, 11, 72]. 

We hope that this paper, in addition to providing a simple theoretical basis for 
treating a new class of nonstationary dispersive media, will demonstrate the utility 
of explicitly recognizing kinetic constraints and of maintaining a clear distinction be- 
tween kinetic and dynamic aspects of wave behavior. Section 5.2 is dedicated to an 
exploration of the interplay between dispersion, nonstationarity, and inhomogeneity 
from a strictly kinetic perspective. We introduce three basic kinetic relationships 
representing three basic preserved symmetries. Along with appropriate dispersion 
relationships, they allow for the derivation of all the subsequent kinetic results. Al- 
though these relationships have not, to our knowledge, been explicitly presented to- 
gether before, they have been implicit in many previously derived results. We show 
that they lead to Snell's Law for a motionless interface and to appropriate Doppler 
shifts for a moving one, that they lead to the preservation of wavelength for homo- 
geneous nonstationary media and to the preservation of frequency in inhomogeneous 
stationary media, how they give rise to Biancalana et a/.'s "generalized frequency" 
and give insight to Avitzour and Shvets' proposal for compressing pulse spectrum 
without altering carrier frequency. We also find some effects that are described here 
for what may be the first time. For example, we find that dispersion modulates the 
frequency response of a wave to temporal changes in the refractive index and that 
Doppler reflections may lead to large changes in pulse bandwidth when the moving 
interface approaches the group velocity. We emphasize that although we have de- 
rived these results with controllably dispersive nonlinear optical media in mind, the 
results apply to any type of wave propagating through any linear medium. 

Unlike kinetic effects, the dynamic aspects of wave behavior depend on the spe- 
cific microscopic details of interactions between nonstationary media and fields. In 
Section 5.3, we develop boundary conditions appropriate to an idealized version of 
the nonlinear optical media currently used to achieve exotic dispersion. Unsurpris- 
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Figure 5.1: An incident plane wave with k-vector k t interacts with a planar interface. 
This results in a refracted wave with k-vector k r and a transmitted wave with k-vector 
k t . For a non-dispersive and isotropic medium, the refractive index for the reflected 
wave (n r ) is the same as that for the incident wave {ni). 



ingly, we find that these boundary conditions are fundamentally different from those 
posited for nondispersive media. We compare pulse transformations wrought by a 
changing phase velocity in a specific, homogeneous, nondispersive medium (originally 
described by Morgenthaler [12]), with those wrought by a change in group velocity 
in our idealized dispersive medium. We note many interesting similarities and dif- 
ferences in their effects on twenty key quantities such as energy, momentum, and 
photon density. 



5.2 Kinetics of dispersive, inhomogeneous, non- 
stationary media 

Consider a plane wave incident on an infinite planar interface, as drawn in Fig. 5.1. Its 
wave vector and frequency describe periodic translational symmetries in 3 dimensions 
of space and 1 of time. A planar interface breaks symmetry in only one dimension, 
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i.e., the one parallel to its normal. If an interface is motionless, then its normal is 
purely spatial; in the geometry of Fig. 5.1, we could write h = x. The normal of 
a moving interface in 3+1 dimensions is partially rotated into time; if we introduce 
the temporal unit vector ct and a temporal rotation angle 0, then we may write 
n = cos0x + sin0ct. The angle <fi is related to the velocity v and the normalized 
velocity f3 by 

- = /3 = tan0. (5.1) 

c 

All the other translational symmetries of the incident plane wave are preserved under 

reflection and transmission. Thus, 

Kyi fiyr ™yti V / 



">zi — Kzr — i^ztt I "-*•") 



and 



— - k xi (3 = — - k xr (3 = — - k xt (3, (5.4) 

c c c 

where the geometry is again defined by Fig. 5.1. These three equations, combined 

with appropriate dispersion relationships, provide a basis for the kinetic aspects of 

wave behavior. 

Equations (5.2), (5.3), and (5.4) represent quantities that are unchanged by a 
planar interface. Extending this concept, we see that if these quantities are not 
altered by a single interface, they cannot be altered by multiple similar interfaces 
("similar" means here that they share the same velocity and the same normal). By 
progressively increasing the number and variety of similar interfaces, we can similarly 
see that the same quantities are preserved under interaction with any pattern of 
medium parameters that obeys the traveling wave law, 

f = f(x-vt). (5.5) 
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When the planar disturbance described by the function / in Eq. (5.5) is more 
complicated than a single interface, the description of the wave response in terms 
of a single incident, reflected, and transmitted wave becomes inadequate. In order 
to express conserved quantities appropriately for this more general case, we now 
introduce the convention that a line appearing over a given quantity means that the 
quantity is invariant under planar changes in the material parameters. Using this 
convention, we extend Equations (5.2), (5.3), and (5.4) to 

k y , (5.6) 



k z , (5.7) 

and 



— cos(</>) — k x sin(c/>). (5.8) 

Despite their lack of an explicit equality, we refer these statements equations in 
acknowledgment of the fact that they are, in fact, compressed equalities. 



5.2.1 Snell's Law and Doppler shifts 

For a motionless interface, (3 = and Eq. (5.4) reduces to the statement that fre- 
quency is preserved: 

uji = u r = ujf (5.9) 

Combining this with Eq. (5.3) and the expansion k z = nusin8/c yields 

sin 8i = sin 8 r = sin 6 t , (5.10) 

which is Snell's Law (note that Snell's Law applies to the reflected wave as well as 
the transmitted wave). 
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If we allow the interface to move at a velocity v , frequency is no longer conserved, 

and Snell's Law will no longer hold 2 . If both sides of the interface are dispersionless, 

then a simple extension of Snell's Law obtained by combining Eq. (5.4) and Eq. (5.3) 

may be useful: 

n t sin 9t Uj sin 9j , . 

1 — nt cos Ot/3 1 — rii cos 9i(3 

One interesting fact about this formula is that it shows that there are some cases 

where there two sets of angles and frequencies for the transmitted wave 3 . When 

dispersion is important, Eq. (5.11) loses its utility because 9 t and n t both become 

frequency dependent. We may still find angle, refractive index, and frequency using 

both relationships 

n t (u t )u t sin 9 t = m(uji)uji sin 9 h (5.12) 

and 

u t (l -n t (uJt) cos 9 t j3) = 0^(1 -n^Wi) cos 9ij3), (5.13) 

combined with an explicit function for n t (oJ t ) to calculate the properties of the trans- 
mitted wave. We will explore the effects of dispersion in a simplified case, that of 
the transmitted wave at normal incidence, in the next subsection. 

The relationships for the reflected wave, 

n r {uj t )u) r sin9 r = n^u^Ui sin 9 iy (5-14) 

and 

u> r (l — n r (u r ) cos9 r (3) = u>i(l — riiiuii) cos 9i(3), (5.15) 

2 For a moving object in a vacuum, this is easily solved via Lorentz transformation. For 
general nonstationary objects, such a transformation may be costly: isotropic media may 
be come anisotropic. In addition, Lorentz transformations cannot be used to simplify the 
treatment of interfaces moving faster than the speed of light. 

3 Multiple transmitted solutions can correspond to the absence of a reflected solution. 
This has been noted in the past for moving interfaces. See [62] 
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combined with an explicitly function for n r (u r ) appear identical to those for the 

transmitted wave. However, the behavior is different in this case because the medium 

that defines n r is the same as the medium that defines n^ and because the sign of 

cos# is different in the reflective case. Numerical investigation suggests that the 

Doppler effect grows large when \v g \ ~ v for an approaching wall. In Subsection 

5.2.5 we will show analytically that this is the case for a wave that is normal to 

the interface. For now, we note that for normal incidence on a moving interface in 

a vacuum, cos8 r = — 1, cos#j = 1, and n r = rii = 1. Then Eq. (5.15) gives the 

standard reflective Doppler shift, 

1-/3 . 

uj r =Ui—— — . (5.16) 

Note that in the geometry of Fig. 5.1 a positive value for means that the interface 
is receding. 

5.2.2 Transmission at normal incidence 

In controllably dispersive media, large changes in group velocity may be accompanied 
by minute changes in phase velocity. In addition, the phase velocity changes are 
gradual in the sense that they occur over several carrier frequency cycles. Under 
these conditions, reflection is negligible. At normal incidence, cos 8i = cos #4 = 1. 
Then Eq. (5.8) simplifies to 



— cos0 — ksm(fi, (5-17) 

c 

where we have removed the subscript x from k x in acknowledgment of the fact that at 

normal incidence k = k x . This mixture of wave vector and frequency is the conserved 

quantity that Biancalana et al. found and referred to as the generalized frequency 

71]. When = 0, the medium is stationary and frequency is preserved. When 

cf) = ±7r/2, the medium is homogeneous and wave vector is preserved. In between 

these extremes, a mixture between wave vector and frequency is preserved. 
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The preservation of the generalized frequency is accompanied by a second con- 
straint: 

^ = v p , (5.18) 

where v p is the phase velocity. When the phase velocity changes, the ratio uj/k 
must change as well. Together, Eqs. (5.17) and (5.18) determine the frequency and 
wave vector responses. We now look at the impact of dispersion upon these responses. 
Consider a medium for which the spectral dependence and the time dependence of the 
refractive index are independent so that we may apply n = n (x, t) + (5n/(5uj))(u — 
uj ) for frequencies near u . From 

(5.19) 



d 

dn 


( — cos cf) — k sin <p j = 


= 0, 


we obtain 








duj 


wsin0 






dn 


cos — rig sin ' 




where n g ■ 


= n + co5n/Soj. 


Since 


tan 


du 


uv/c 







v/c, we can also write 



(5.20) 



(5.21) 



dn 1 — v/vg ' 

which emphasizes the large frequency shifts expected when the interface velocity v 
approaches the group velocity v g . 

5.2.3 Preservation of a generalized bandwidth in 1+1 di- 
mensions 

The group velocity is defined using v g = du/dk. When group velocity dispersion 

may be neglected, either because it is small, or because we are limiting ourselves to 

a narrow spectrum, the approximation 

Au 

U « v. (5.22) 
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becomes useful. As formalized in Eq. (5.22), when a narrow band pulse experiences 
a change in group velocity, either its spectral bandwidth Aw, or its spatial bandwidth 
Ak, or both, must change. Extending Eq. (5.17), we find a generalized bandwidth, 



cos <p — Ak sin 0, (5.23) 

c 

that is preserved when a pulse interacts with an interface moving at the velocity 
defined by and Eq (5.1). Eq. (5.23) shows that the division of change between 
spectral and spatial bandwidth depends on their comparative magnitudes and on the 
interface velocity. 

The preserved generalized bandwidth provides us with a simple tool that we can 
use to analyze one particular aspect of the stopped light experiment of Liu, Dutton, 
Behroozi, and Hau [ ]. In the Liu experiment, a coupling beam was used to control 
the group index perceived by a probe beam in real time. However, this coupling 
beam had a finite speed and was copropagating with the probe beam. Because of 
this, changes to the group index perceived by the probe beam were not instantaneous, 
but moved at nearly the speed of light. In their report [ ], and in earlier theoretical 
work [1 3], the speed of the coupling beam was treated as effectively infinite, because 
it was so much faster than the slowed advancement of the probe pulse. Using the 
preserved generalized bandwidth, we now show that treating the velocity of the 
coupling beam as essentially infinite was justified. 

When the pulse described by Liu et al. entered the slow light medium, it crossed 
a spatial interface, and <fi was 0. The preserved quantity associated with this case 



was Aoj/c, meaning that spectral bandwidth was preserved. The spatial bandwidth 
(Ak) and the group index (n g ) were increased by 7 orders of magnitude, and Eq. 
(5.22) remained satisfied. 

While the new foreshortened pulse was propagating through the slow-light medium, 
the coupling beam was then turned off. The speed of the coupling beam was ap- 
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proximately c, so that <fi was approximately 7r/4. Thus, the preserved generalized 



frequency was approximately (Aw/c — Ak)/y2, which appears to distribute changes 
evenly between spectral and spatial bandwidths. However, the slow speed of the 
probe envelope was reflected by the fact that Ak was larger than Aw/c by a factor 
of 10 7 . According to the preservation of the generalized bandwidth, Au/c would 
have gone to zero while Ak changed by 1 part in 10 7 . This is the size of the error 
introduced by assuming that the coupling speed was effectively infinite. 

The end result of the Liu experiment was a reconstitution of the original pulse. 
In the next section we use the concept of preserved generalized bandwidth to show 
how the pulse may be spectrally compressed or expanded using transmission in con- 
trollably dispersive media. 

5.2.4 Pulse compression/decompression using controlled dis- 
persion 

The generalized bandwidth is conserved across any number of interfaces so long as 
they have the same velocity. By mixing interfaces of different velocities, we can 
change that bandwidth. In Fig. 5.2, a pulse begins in a medium with a group index 
n g of 1, transitions to one with n g = 2, and then back to n g = 1. When both 
transitions are spatial (S corresponds to a motionless planar interface) or when both 
are temporal (T corresponds to time-dependent but spatially independent change to 
the medium-the infinite velocity case) the pulse ends up in its original form. By 
mixing spatial and temporal transformations, the pulse can either be compressed 
or extended. In this way the pulse bandwidth may be compressed or decompressed 
while the pulse envelope is scaled longitudinally but otherwise unchanged. 

This mechanism for pulse compression and decompression was recently proposed 
and demonstrated in simulation for a plasma geometries, where the group velocity 



71 



Chapter 5. Nonstationary electromagnetics of controllably dispersive media 




Figure 5.2: Field energy density verses spatial extent of a pulse. S: when a pulse 
undergoes a spatial change in group index (n 9 ) it scales longitudinally but not tem- 
porally. T: when a pulse undergoes a temporal change in (n 9 ), it changes temporally 
but not longitudinally. By mixing spatial and temporal interfaces, it is possible to 
scale the pulse bandwidth and longitudinal extent without otherwise changing the 
envelope shape. 



may be controlled dynamically via an applied magnetic field [ ]. We note here 
that the Liu experiment demonstrates all the capabilities necessary to perform pulse 
compression or expansion using nonlinear optics [9]. 

5.2.5 Dispersive modulation of the 1-D Doppler shift 



In a dispersive medium, Doppler shifts are difficult to predict intuitively because of 
the interdependence between refractive index and frequency. Here we find an intuitive 
formula for the change in frequency in the case where the group velocity differs from 
the phase velocity, but where the range of frequencies concerning us is sufficiently 
narrow that we may neglect group velocity dispersion. Expanding Eq. (5.4) for the 
case of reflection at normal incidence, we get 



UJ; 



71; 



U r 



1 + n r 



(5.24) 
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To derive a Doppler shift, we must introduce an explicit dispersion relationship 
n = n(u>). In general, the resulting equation is analytically intractable though a 
numerical solution is readily found. An interesting argument for the physical rele- 
vance of the concept of group velocity to Doppler shifts in dispersive media is that 
if we require that the group velocity be constant with frequency a simple analytical 
solution results. Requiring that n g , the group index, remain constant with frequency 
gives a refractive index of the form n = n g + b/u, where b can be any real constant. 
Applying this dispersion relationship to Eq. (5.24) yields the exact expression, 

A v , 

- = -2n-, 5.25 

ui c 

where (3 = v/c, A = u r — u>i, u r + u>i = 2cu, and n = n{uj). Eq. (5.25) is attractive 
because of its symmetry but incomplete because it depends on the quantity u, which 
it does not resolve. Rewriting Eq. (5.25) in terms of Ui and A 

* = - 2 " (C ^ /C , (5.26) 

Ui 1 + V/Vg 

which clearly reveals the role of dispersion in determining the expected Doppler shift. 
Eq. (5.26) has the interesting interpretation that the Doppler shift will be large when 
the interface is approaching (if v g > 0) or receding (if v g < 0) with a speed \v\ ~ \v g \. 
This result is exact when the group velocity is constant with frequency. 

5.2.6 Dispersion and the 1-D Doppler effect on pulse band- 
width 

Intuitively, we expect that a pulse reflecting from an interface whose velocity is close 
to that of the pulse will be temporally stretched if the interface is receding and 
compressed if the interface is approaching. We expect this temporal stretching/com- 
pressing should have a corollary decrease or increase of the spectral bandwidth, Aw. 
We now show that this is the case. 
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Adapting expression (5.23), we find 

Auji - Ak { v = Aw r + Ak r v, (5.27) 

where Aw», Aki, and Au r , Ak r are the spectral and spatial bandwidths of the incident 
and reflected pulses, and v is the interface velocity. Assuming that the initial and 
final pulses are sufficiently narrow such that each may be associated with a specific 
group velocity, we may combine cv 9 = Au/Ak with Eq. (5.27) to obtain 



u gr 



where v gi is the group velocity of the incident pulse and v gr is the group velocity of 
the receding pulse. As always, a positive sign for v means a receding interface and 
a positive sign for v g means a positive group velocity. Thus, as expected, when a 
pulse reflects from an interface that is receding at roughly the pulse group velocity 
of the approaching pulse, the reflected pulse has a narrowed bandwidth. When a 
pulse reflects from a boundary that approaches at roughly the group velocity of 
the reflected pulse, the reflected pulse is temporally compressed and has a broad 
bandwidth. 

5.2.7 How group velocity modulates the frequency response 
of a pulse to a time-dependent refractive index 

When interface velocity, v , is zero, <fi is also zero and u is conserved. When velocity 
is infinite, = ±7r/2. The refractive index is a function of time only and k is 
conserved. In this way, the interface model is easily extended to treat the case 
of a time-dependent (but spatially homogeneous) medium. Since k is conserved, 
frequency must compensate for a change in refractive index according to 

Ti 

u(t) = -^u , (5.29) 

n(t) 
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where we take the medium to have an initial refractive index n and the pulse to 
have an initial central frequency u . If n were driven to 0, u would be up-converted 
by many orders of magnitude. The changing refractive index associated with plasma 
generation has been proposed for use in frequency up-conversion in several works 
74, 75, 63, 64, 65, 67, 68]. On a more pedestrian level, frequency response to a 
time-dependent refractive index is also at the heart of the function of acousto-optic 
and electro-optic modulators. 

Interestingly, the frequency response to a changing refractive index can be mod- 
ulated by the group velocity of a pulse in a medium. To see this, we begin with 
a dispersive medium with a nondispersive time-dependence whose refractive index 
takes the form n = n Q (u) + tSn/8t. (This might be taken to model, for example, 
a gas that is a mixture of two atoms, one with resonance that is near u, the time- 
dependent frequency of the wave traveling through the medium, and another whose 
resonances are all spectrally distant but whose concentration changes with time.) 
Then, from Eq. (5.29) we find 

duj u 5n . 

— = — — , 5.30 

dt n g 5t v ' 

where n g is the group refractive index, defined by n g = n + uSn/Sco. One interpre- 
tation of Eq. (5.30) is that the frequency response of a pulse to a temporal change 
in refractive index is proportional to the distance the pulse travels while the change 
occurs. 



5.3 A comparison of two nonstationary homoge- 
neous media: kinetics and dynamics 

In this section we seek to augment our understanding of the effects of a time- 
dependent group velocity in a dispersive medium by comparing them with the effects 
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Number 


Quantity 


Symbol 


Morgenthaler 


Liu 


1 


Phase velocity 


v p 


Vpf/V pi 


1 


2 


Group velocity 


V 9 


Vpf/Vpi 


V 9fl V 9i 


3 


Permittivity 


e 


e /A* 


1 


4 


Permeability 


n 


\ifl\ii 


1 


5 


Wave number 


k 


1 


1 


6 


Spatial bandwidth 


Ak 


1 


1 


7 


Central frequency 


UJQ 


Vpf/Vpi 


1 


8 


Spectral bandwidth 


Auj 


Vpf/Vpi 


V 9fl V 9i 


9 


Total energy density 
Electric displacement 


D 


Vpf/Vpi 

1 


1 


10 


V V 9f/ V 9i 


11 


Magnetic induction 


B 


1 


V V 9f/ V 9i 


12 


Electric field 


E 


e */ e / 


V V 9f/ V 9i 


13 


Magnetizing field 


H 


Hi/fif 


V V 9f/ V 9* 


14 


Photon energy 


Huj 


Vpf/Vpi 


1 


15 


Photon number density 


N 


1 


1 


16 


Poynting vector 


S 


(Vpf/Vpi) 2 


Vgf/Vgi 


17 


Abraham momentum density 


E x H/c 2 


iVpf/Vpif 


V 9f/ V 9i 


18 


Minkowski momentum density 


Dx B 


1 


V 9f/ V 9i 


19 


Canonical momentum density 


N(nhu/c) 


1 


1 


20 


Angular momentum density 


N(h) 


1 


1 



Table 5.1: A comparison of narrow-pulse transformations by two dynamic media. 

of a time-dependent phase velocity in a nondispersive medium. To do so, we now 
define two scenarios, which we call Scenarios A and B. In both scenarios, we imagine 
a spectrally narrow pulse traveling through a time-dependent medium and ask how 
that pulse is transformed as the medium changes with time. 

In Scenario A, we take the medium to be nondispersive, spatially homogeneous, 
and isotropic. We assume that changes in the medium are adiabatic, spatially ho- 
mogeneous, and isotropic. We allow e and \x to vary gradually and independently 
over time. We assume that loss at all times is negligible, meaning that e and \x are 
taken to be real for the narrow spectral range of the pulse. Because the medium is 
nondispersive, the group velocity is equal to the phase velocity. 
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In Scenario B, we take the medium to be dispersive, spatially homogeneous, and 
isotropic. We assume that changes in the medium are adiabatic, spatially homo- 
geneous, and isotropic. We fix e and fi for the central frequency, but allow their 
slopes to vary gradually and independently over time. We take loss at all times to be 
negligible, meaning that e and \x are taken to be real for the narrow spectral range 
of the pulse. 

5.3.1 Kinetics: basic symmetries 

The spatial homogeneity of Scenarios A and B implies that 

k f = ±k h (5.31) 

and the adiabaticity of the two scenarios eliminates reflection, requiring that we use 
only the positive sign. Spatial homogeneity also implies that 

Ak f = Aki (5.32) 

for the pulses of the two scenarios. Since \k\ = nu/c and Ak w n g Au/c, changes 
in the phase and group velocity imply proportional changes in the central frequency 
and pulse bandwidth respectively. Since n g = n + u5n/5co, changes in n imply 
changes in n g for a nondispersive medium and therefore changes in n lead to changes 
in bandwidth as well. 

To continue our comparison, we now specify particular boundary conditions for 
the two changing media. In doing so, we lose the generality of Section 5.2, but gain 
the ability to examine the effects of changes to specific, time-dependent media on 
basic quantities like energy and momentum densities and reflection and refraction 
coefficients. 

We take our two media from previous works. The non-dispersive medium is taken 
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from Morgenthaler [12]. The dispersive medium is an idealization of the dynamic 
EIT medium used in the Liu experiment [9]. 

5.3.2 Dynamics: boundary conditions and reflection 

Morgenthaler's nondispersive medium 

Morgenthaler allowed for the electric permittivity (e) and magnetic permeability (y) 
of the medium to be functions of time, but not of space or direction or frequency. 
That is, e(r, k, u>,t) = e(t), and /z(r, k, cu,t) = //(£). He assumed that charge and 
flux, and therefore the electric displacement (D) and the magnetic induction (B), 
would be conserved at an interface boundary. Under this assumption, the reflection 
and transmission coefficients for the electric field are given by [12] 



E r L (e % [JL& , 



E{ 2 \€ r y /i r e 
and 



Et| L ^ + /M). (5.34) 



|Ej| 2 \€ r y /i r e 
In any Morgenthaler medium, the magnitude of the magnetic field is given by |H| = 
|E|/?7, where n = y/i/e is the impedance. In general, there is a reflected wave. 

We consider the case where changes to the medium are adiabatic and homoge- 
neous. In this case the reflected power vanishes. Using the fact that what is preserved 
across one interface must be conserved across many parallel interfaces (see the be- 
ginning of Section 5.2), we find that we may take the boundary conditions across an 
adiabatic change to be 

D/ = D; (5.35) 

and 

B f = B z . (5.36) 
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A = 



Figure 5.3: EIT energy diagram. The coupling beam and the medium may be taken 
together to form an effective medium as felt by the probe beam. The properties of 
the effective medium depend on the state of the coupling beam. 



The values of E and H at any time can then be found via e and \i. 



Liu's dispersive medium 



The Liu experiment involves a pulse interacting with an EIT medium across a narrow 
spectral range around the EIT resonance. The system as a whole comprises two 
beams of light interacting with a cloud of cold atoms. The two beams are tuned to 
two connected atomic transitions which share an excited state. The stronger beam, 
referred to as the coupling beam in Figure 5.3, is considered in combination with 
the cloud of atoms to form the EIT medium, as seen by the weaker beam, which 
we call the probe beam. When EIT is established, the atoms are in a coherent dark 
state, so-called because it does not couple to the probe beam. The dark state may 
be represented as a combination of electronic levels |1) and |2) as 



\D) 



fi c |l) — fi p |2) exp[i(k p — k c 



r — i (u r 



u c )t] 



(5.37) 



If the probe beam is sufficiently weak such that the atomic number density remains 
significantly larger than the probe photon number density, then the EIT medium 
will be approximately linear in probe intensity, so that the only sizable nonlinear 
effects of the EIT medium are the effects of the coupling beam on the probe beam. 
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Surprisingly, this linear regime may be maintained even as the coupling beam is 
completely attenuated [73]. 

For our nondispersive medium, we idealize the Liu experiment ['. '] by neglecting 
decoherence and other losses and require that total pulse energy be conserved under 
changes to the group velocity. Thus, 

u t f = u t i, (5.38) 

where Ut represents the total energy density associated with the wave including not 

only that portion that is in field form but also that portion that is coherently stored in 

the medium and is a cycle-averaged quantity. For a narrow-band pulse propagating 

through a lossless but dispersive medium, we may define a non-dispersive energy 

density (u n ) which is related to the dispersive energy {u t ) by 

n 
u n = fh—. (5.39) 

n g 

In terms of the macroscopic fields, u n is given, for the isotropic medium we are 
considering, by 

u n = l -{eE 2 + u.H 2 ), (5.40) 

where E 2 and H 2 are time averaged quantities. (Note that u n matches the standard 
form for the electromagnetic energy density in a lossless, dispersionless, macroscopic 
medium [32]. When dispersion becomes negligible, n g « n, and u t reduces to this 
form. ) 

If, as we have assumed for our comparison, changes in the group velocity occur 
around a constant refractive index, then we may take the quantity u n n g to be fixed. 
We will also assume that the impedance at line center is unchanged, so that the ratio 
between H and E remains fixed. Because the refractive index at line center will also 
remain fixed, all macroscopic electromagnetic fields scale together. Therefore, all 
terms quadratic in any pair of fields scale with the group velocity. The direction of 
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the fields is unchanged by the homogeneous change in the isotropic medium. We 
summarize this as 

F/v^/ = Fi^V, (5.41) 

where F could be any of the four macroscopic electromagnetic fields. 

5.3.3 Transformations wrought under changing phase and 
group velocities in the 'Morgenthaler' and 'Liu' media 

We are now able to perform a more complete comparison between scenarios A and B. 
In Table 5.1, we list 20 different quantities and show how they are scaled as the media 
of Scenarios A (Morgenthaler) and B (Liu) change. Quantities 1-4 summarize the 
difference between the two scenarios, emphasizing that the phase velocity and group 
velocity change in Scenario A while only the group velocity changes in Scenario 
B. The scaling of quantities 5-8 depends only on symmetry (see Section 5.2) and 
therefore are general consequences of the scaling of quantities 1-4. The scaling of 
quantities 9-13 may be taken as expressions of the boundary conditions for the two 
media. 

The remaining quantities (14-20) can be found in simple ways from a knowledge 
of the scaling of the first 13 quantities. The per photon energy (14) scales with 
frequency, which scales as in line 7. That the photon number density (15) is con- 
served is related to our assumption of no loss. The Poynting vector (quantity 16) 
scales as the square of the phase velocity for the non-dispersive case. This quadratic 
dependence combines the linear dependence of the total energy density on the phase 
velocity with the fact that the group velocity is equal to the phase velocity for a 
nondispersive medium. Because the dispersive case assumes that the phase velocity 
remains unchanged, the Poynting vector of the Liu medium feels only the effect of 
the change in group velocity. The Minkowski momentum (18) is conserved in the 
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nondispersive case, but is proportional to the group velocity for the nondispersive 
case. This second dependence underscores the fact that the Minkowski momentum, 
as defined by D x B, does not give the total momentum when dispersion is taken 
into account. The canonical momentum (19), which does give the total momentum, 
is conserved for both media. Because the total photon number is conserved for both 
media, so is the total angular momentum density associated with the pulse for both 
media. We note that although the total angular momentum does not change, its 
distribution between field and medium does. 

In concluding this section we emphasize three points. First, to understand the 
effects of a particular medium on a pulse, it is necessary to specify appropriate bound- 
ary conditions that depend on the details of the way in which the changing medium 
interacts with the field. Second, boundary conditions that are reasonable for a non- 
dispersive medium (Scenario A) may fail for a simple dispersive medium (Scenario 
B) for the particular reason that changes in dispersion imply a shifting of energy 
from field to medium. It then becomes necessary to model different changes on a 
case-by-case basis. Finally, a general comparison of the scaling of different quantities 
for Scenario A and Scenario B shows that the changes wrought in a controllably 
dispersive nonstationary medium are qualitatively different from those wrought in 
a nonstationary nondispersive medium. In this sense, controllably dispersive media 
open up a new and potentially fruitful niche in nonstationary electromagnetics. 



5.4 Conclusion 

Controllably dispersive media are an experimental reality [8, 9, 57]. All the require- 
ments necessary to perform interesting nonstationary experiments in controllably 
dispersive media have been demonstrated for one experimental system [9, 57], are 
near at hand for another [8], and appear to be feasible for others [15, 11]. In this 
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paper we have outlined a basic theory for the treatment of controllably dispersive 
media from the perspective of nonstationary electromagnetics. 

Although allowing the properties of a medium to depend not only on position 
but also on both time and frequency opens up a large parameter space, the behavior 
of waves in this space is subject to basic kinetic constraints imposed by symmetry. 
We have explicitly defined the requirements of symmetry for the case of a plane wave 
interactive with a moving planar interface. We have justified these requirements by 
rederiving well-established results, such as Snell's Law and the free space Doppler 
shift, and used them to simply derive analytical descriptions of more esoteric phe- 
nomena, such as the effect of dispersion on the reflective Doppler shift. Using the 
outlined kinetic theory, we have provided an analysis that corroborates a recent pro- 
posal [1 1] for the compression and decompression of pulse envelopes in magnetized 
plasmas using strictly nonstationary effects and have shown that the proposed mech- 
anism may be applied in any controllably dispersive medium. We have introduced 
a simple expression that shows how a pulse may be significantly compressed or de- 
compressed under Doppler reflection-even if the Doppler shifts themselves are small. 
Finally, we have shown that dispersion modulates the frequency response of a wave 
to a temporal change in the refractive index-the response is damped for slow light 
media, and amplified for fast light media. 

Unlike kinetic effects, which rely only on simple symmetries, dynamic effects 
are determined by boundary conditions which depend upon particular details of the 
medium at hand. We have compared dynamic effects in two simple, homogeneous, 
adiabatically nonstationary media-one which is nonstationary in phase velocity and 
a second which is nonstationary in group velocity. We compared the effective bound- 
ary conditions associated with these media and the changes wrought in spectrally 
narrow pulses propagating through them. We have found substantial differences in 
both the boundary conditions and the effect of the nonstationary medium on funda- 
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mental quantities such as energy, and electromagnetic momentum. These differences 
suggest that it may be useful to consider controllably dispersive media as a new and 
potentially fruitful category of nonstationary electromagnetic media. 
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Chapter 6 

Dispersion and the quantum noise 
associated with a single cavity 

mode 



Abstract 

Anomalously dispersive cavities, particularly white light cavities, have been proposed 
for use in LIGO-like gravity wave detectors and in ring-laser gyroscopes. In this paper 
we analyze the quantum noise associated with anomalously dispersive cavity modes. 
The vacuum field energy associated with a particular cavity mode is proportional to 
the cavity-averaged group velocity of that mode. For anomalously dispersive cavities 
with group index values between 1 and 0, this means that the total vacuum field 
energy associated with a particular cavity mode must exceed Huo/2. For white light 
cavities in particular, the group index approaches zero and the vacuum field energy 
of a particular spatial mode may be significantly enhanced. We predict enhanced 
spontaneous emission rates into anomalously dispersive cavity modes and broadened 
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laser linewidths when the linewidth of intracavity emitters is broader than the cavity 
linewidth. 



6.1 Introduction 

When the finesse of an optical cavity is improved, buildup increases but resonance 
bandwidth decreases. Resonance bandwidth may be increased without degrading 
buildup by decreasing cavity length. For some applications, such as interferometry- 
based gravity-wave detection and Sagnac interferometry, decreasing cavity length 
degrades overall system performance. For other applications, such as fast single- 
photon generation [76, 17], decreasing cavity length may be profitable, but is not 
possible past a minimum length. 

When actual decreases in length are either unprofitable or impossible, it is natural 
to consider the use of anomalous intra-cavity dispersion. For a given cavity geometry, 
buildup scales with finesse (F) 1 , while resonance bandwidth scales inversely with 
quality (Q). Quality and finesse are related by Q = F * (y L g /c), where v is the 
resonant frequency, c is the speed of light, and L g is the round-trip group optical 
path difference. L g itself is defined in terms of the round trip phase <fi and the angular 
frequency u by L g = cdtfi/du [ ]. In the presence of dispersion, L g = £n g , where £ is 
the cavity round trip length and n g is an effective cavity group index. The response 
of the ratio F/Q to changes in n g is identical to its response to changes in £. Unlike £, 
however, n g may approach zero for a particular range of frequencies [78, 31]. Cavities 
where n g approaches zero at a cavity resonance have been referred to as "white-light 
cavities," in reference to their dispersion-broadened resonances [79]. 



x By buildup we refer to the resonant enhancement of intracavity intensity as normalized 
by the incident intensity of a mode. For high- finesse cavities where the round trip electric 
field transformation is E — > Ee~ l ^{\ — 5), the buildup is A/5 2 . 
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The white-light cavity literature is centered around two potential applications, 
gravity wave sensing [ 7 9, 80, 81] and rotation sensing [I ]. LIGO-like gravity wave 
detectors take the form of large Michelson interferometers. In signal recycling, mir- 
rors may be placed near the entrance to each Michelson arm, forming a pair of 
Gire-Tournois cavities. In power recycling, a mirror is placed in advance of the cen- 
tral beam splitter, making the entire interferometer into an optical cavity. In either 
case, the increase in signal due to enhanced buildup is gained only at the expense of 
bandwidth, and white-light cavities have been considered as a way to combine large 
bandwidth with large buildup. 

Rotation in a ring-resonator gyroscope is equivalent to cavity elongation/trunca- 
tion when resonator and mode are co-/counter-rotating [82]. The frequency response 
of a cavity mode to a change in length may be enhanced via anomalous dispersion; 
as the round-trip length £ changes, the resonant frequency u associated with a 
given mode changes according to du /d£ = —nu /(£n g ), which grows large as n g 
approaches 0. Shahriar et al. found that the increase of frequency sensitivity in 
a passive anomalously dispersive ring cavity would be counterbalanced by the con- 
current increase in resonance bandwidth, leaving no net gain in rotation resolution 
[82]. However, they argued that this cancellation could be avoided by using an active 
interferometer in the form of a white-light ring-laser gyroscope. 

The effectiveness of white-light cavities for any sensing application is partially 
determined by white-light cavity noise. Some implement at ion- dependent sources of 
noise have been quantified by Wicht et al. [ i] for a double-lambda system and by 
Sun et al. \ i] for a double gain system. A comparison of several implementations 
can be found in a later work by Wicht et al. [ ]. However, the white-light cavity 
literature does not include, to our knowledge, any work on the noise due to vacuum 
energy intrinsic to a white-light cavity. 

In this paper, we show by simple and largely classical arguments that the vac- 
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uum field energy associated with the ground state of an anomalously dispersive cavity 
mode diverges as that mode approaches the ideal (n 9 — > 0, F — > oo) associated with 
a white-light cavity. This divergence is intrinsic to the definition of a white-light 
cavity and is independent of any specific implementation. It has the consequence 
that spontaneous emission of an excited particle into the white-light cavity mode 
may increase substantially. This, in turn, has two consequences. First, the width 
of a white-light laser line must be broadened by anomalous dispersion. This broad- 
ening, if applied to the configuration proposed by Shahriar et al. would cancel the 
increased frequency sensitivity of a white-light ring-laser gyroscope. Second, it sug- 
gests that anomalous dispersion may be useful in enhancing the quantum yield of 
single photon emitters where radiative decay into the mode of interest must compete 
with nonradiative decay or with radiation into other spatial modes. 

In the following section we review the result of applying the standard classical ex- 
pression for the narrow-bandwidth energy density in a lossless, dispersive medium to 
the quantization of a closed cavity. We introduce a distinction between two energies, 
which we call the "field energy" and the "total energy," and which are related by the 
ratio between cavity-averaged phase and group velocities. The total cavity energy 
is quantized as a simple harmonic oscillator, but only the field energy interacts with 
dipoles inside the cavity. In Section III, we explore the vacuum field noise in an 
anomalously dispersive cavity using an alternative coupled-cavity approach. This 
approach allows us to find an expression for the vacuum fields that is independent of 
any specific expression for the electromagnetic energy density in a dispersive medium. 
It also provides a simple view of the effect of dispersion on the power spectrum of 
the vacuum field energy. We show that if we apply the assumptions of Section 6.2 
(no loss and negligible group velocity dispersion for the spectrum of the resonance), 
then the results of Section 6.2 follow. We then relax these assumptions using a nu- 
merical model of a physical medium and show numerically that as the round trip loss 
approaches zero, the quantum field noise of a white light cavity resonance diverges. 
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This occurs even if the direct effect of the gain of the anomalously dispersive medium 
is neglected. We briefly address the relationship between this increased noise and 
the quantum limited laser linewidth of anomalously dispersive cavities. 



6.2 Spontaneous emission in a lossless, dispersive 
medium 

Electromagnetic field quantization in an evacuated, closed and lossless cavity pro- 
ceeds as follows: 

1. A complete set of orthonormal electromagnetic modes is identified. 

2. The volume-integrated energy of each of these modes is represented. 

3. Each mode is quantized as though it were a simple harmonic oscillator. 2 

When a dispersive medium is introduced into a cavity, both steps 1 and 2 become 
problematic. Step 1 becomes problematic because a dispersive medium must have 
loss (or gain), which couples the electromagnetic fields of a dispersive cavity to 
external degrees of freedom. True electromagnetic modes can then only be found by 
explicitly including these external degrees of freedom. 

Step 2 becomes problematic because dispersion also complicates the concept of 
electromagnetic field energy density. Whether or not energy that has been stored 
by a dispersive medium returns to field form depends on future interactions between 



2 Note that this is the only step that introduces a non-classical element. It is the two 
classical steps that define the relationship between mode excitations and electromagnetic 
fields. 
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medium and field 3 . In other words, the electromagnetic energy density of a dispersive 
medium is nonlocal in time. 

These difficulties with steps 1 and 2 can both be avoided if, rather than seek 
a complete set of orthonormal electromagnetic modes we focus on a single cavity 
mode. Although a causal dispersive medium must have loss (or gain), causality does 
not prevent this loss from being small or even vanishing for particular frequencies. 
Thus, although a causal dispersive cavity must be open in general, it may be closed 
at particular resonances. In addition, a perfect, closed cavity mode has a vanishing 
spectral width. Although the electromagnetic energy density in a dispersive medium 
is, in general, an ill-defined quantity, it is well defined in the absence of loss for 
spectrally narrow excitations. 

6.2.1 Energy density of a quasimono chromatic planar wave 
in a lossless, dispersive medium 

The effect of dispersion on the electromagnetic energy density of a quasimonochro- 
matic excitation in a lossless linear medium is particularly simple for planar waves. 
When only radiative energy transportation is non-negligible, Poynting's theorem 
takes the form 

V ■ S = -du/dt, (6.1) 

where S is the Poynting vector and u is the electromagnetic energy density. If absorp- 
tion, scattering, and group velocity dispersion are all negligible, a direct calculation 
of the divergence of a spectrally narrow planar wave propagating through a dispersive 
medium yields 

o( ,, t) = M, (6.2) 



3 a concrete application of this concept may be found in [84] 
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where v 9 is the group velocity of the planar wave, and the angle brackets denote a 
cycle average. 

If dispersion were neglected, (S)/v 9 in Eq. (6.2) would become (S)/v p , where v p 
is the phase velocity. In other words, if Uf is an expression for the nondispersive 
energy density, then u can be found from Uf via 

u = —Uf, (6.3) 

n 

where n g , and n are the refractive index and the group index. The standard form 
for the cycle averaged electromagnetic energy density in a lossless, isotropic, linear, 
dispersionless medium is [85, 32] 

u, = ^ + *£>. (6.4) 

Substituting Eq. (6.4) into Eq. (6.3) gives 

For spectrally narrow excitations in lossless media, dispersion has a simple rela- 
tionship to energy storage by the dispersive medium. If we call Uf the field energy 
density and u the total energy density then we can also define a third quantity, u s 
as the stored energy density, or the difference between the total energy density and 
the field energy density. In normally dispersive media, u > Uf, and u s is positive 
because electromagnetic energy is dispersively stored in the medium. In anomalously 
dispersive cavities, where < n g < n, u < Uf and u s is negative. Just as u s > 
in a dispersive medium signifies temporary absorption, u s < in an anomalously 
dispersive medium signifies temporary emission. In either case, energy exchange is 
governed by the interaction between the instantaneous spectrum of the planar wave 
and the response function complex susceptibility of the dispersive medium [10]. In 
the limit of a white light cavity, u s = —Uf and u — 0, meaning that all of the field 
energy has been donated by the medium. 
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The energy density given by Eq. (6.2) has the strength that the effect of dispersion 
on the energy density is explicit and simple. However, it was derived for planar waves 
and does not apply when there is interference between waves propagating in different 
directions. 

6.2.2 Energy of a quasimonochromatic mode in a lossless, 
dispersive medium 

A more general but less intuitive form for the energy density may be found by using 
vector identities and Maxwell's equations to make the transformation 

v.b-.» + h.« (,a) 

and then using the narrow-band nature of the excitation and linearity to obtain 

[85, 32] 



u = Re 



d{ue) 
du 



(E 2 ) 



+ Re 



d(u)n) 
du 



(H 2 ) 



(6.7) 



where the angle brackets denote a cycle average. 



Eq. (6.7) shows that the simple relationship between dispersion and energy den- 
sity given for planar waves in Eq. (6.3) does not generalize to non-planar waves. For 
example, consider the energy density of counter-propagating waves of equal ampli- 
tude and frequency and identical polarization in a dielectric. If we assume that e 
and \x and their first frequency derivatives are real, then we can use the identities 
n = ^Jeu./{e u. ) and n g = n + udn/du to obtain 

n„ 1 dine 1 din a , 

y_ _ i I _i_ i_ la q\ 

n 2<ilna;2(ilna; 

Since the medium is a lossless dielectric, ude/du = rje, and udfi/du = 0, where 
n is a real number. Applying these values to Eq. (A. 19) gives n g /n = 1 + n/2. 
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In the standing wave that results from the counter-propagating waves, electric field 
nodes correspond with magnetic field antinodes. At electric field nodes, E 2 = 0, and 
u/uf = 1. At electric field anti-nodes, H 2 = 0, and u/v,f = 1 + rj. Thus, Eq. (6.4) 
does not generalize to non-planar waves. 

However, an analogue to Eq. (6.4) does apply to cavity modes. Assuming loss- 
lessness and vanishing imaginary derivatives for e and fi, we rewrite Eq. (6.7) as 



We now note that the electric and magnetic contributions to the total modal 
energy of a nondispersive cavity mode are equal. That is, for the field associated 
with a mode that satisfies closed boundary conditions and the Helmholtz equation, 

^m d , r ^f^m d ,,, (6 . 10) 



This is equivalent to the fact that the cycle-integrated momentum and position en- 
ergies of a simple harmonic oscillator must be equal. The spatial field distribution 
of a mode, as governed by the Helmholtz equation and a particular set of boundary 
conditions, does not change when dispersion is added, so long as e and fi are not 
changed for the modal frequency. Thus Eq. (6.10) remains true when dispersion is 
taken into account (although it no longer equates the total electric and magnetic 
contributions to the modal energy). 

Applying Eq. (6.10) to an integral of Eq. (6.9) over the modal volume of a ho- 
mogeneous cavity gives 

( 1 dine ldlnu\ f ,, . 

U=[l + -— + —rr^ ~ ) u f d 3 r, 6.11 

V 2d\noo 2d\nou) J } y ' 

where U is the total modal energy and Uf is the field energy density given by Eq. (6.4). 
Because we are taking \x and e and their first derivatives to be real, we can substitute 
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Eq. (A. 19) into Eq. (6.11), which gives 

Tl 

U=-*-Uf, (6.12) 

where the field energy Uf is defined by 

U,= im + ?m*r. (6.13) 



2 2 

Eq. (6.12) is a modal energy analog to Eq. (6.3) for planar waves. 

Uf is a cycle-averaged expression for the energy of a nondispersive electromagnetic 
mode. Eq. (6.12) gives a simple description of the effect of dispersion on the total 
modal energy associated with a given field strength. 



6.2.3 The quantized electric field 

Assuming that a dispersive cavity mode may be quantized as a simple harmonic 
oscillator, we write its Hamiltonian H in terms of an annihilation operator (a) and 
a creation operator (a*) as 

H = hw(a j a+-j . (6.14) 

We note here that the energy associated with H cannot be negative and cannot 
therefore be related to Eq. (6.12) when the ratio n g /n is negative 4 . The electric field 
may then be written in terms of a and a) in the form 

E(r, t) = £af(r)e iuJt + £*a?f* {r)e~ iul \ (6.15) 



4 Interestingly, a second argument suggests the impossibility of a lossless resonance at a 
frequency where there is a negative group velocity. Causality dictates that the combination 
of losslessness and a negative group velocity is only seen in an active medium. If there is 
a frequency with a negative group velocity, then nearby are at least two frequencies that 
share the same wavelength as the on-resonance frequency and that are therefore also at 
resonance. However, it appears that at least one of these must also be at a frequency where 
there is gain. Net round trip gain at resonance in a cavity is not tenable in the steady 
state: the cavity must either lase or have its quality destroyed. 
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where f (r) gives the spatial distribution and polarization of the electric field and 8 
is a complex number whose amplitude gives a characteristic electric field strength. 
To make the effect of the mode volume (V) on the electric field explicit, we choose 
to normalize f according to the rule j f 2 d 3 r = V. Combining Eqs. (6.15), (6.14), 
(6.13), (6.12), and (6.10) then gives 



\£\ = \— Vr 6.16 

" erigV 



hw 



If we restrict our consideration to dielectrics, then e = n 2 , and \£\ — . , 
which is the expression for the electric field strength given by Garrison and Chiao 
[86] who follow Milonni [1 i]. This expression agrees with a more general and earlier 
one implicit in work by Drummond [13]. In addition, dispersive expressions for 
characteristic field strengths (electrical or otherwise) may be related to expressions 
that neglect dispersion [87] using the rule \S diS p ersive \/\S nondispersive \ = y/n/n g . 

Eq. (6.16) suggests three separate ways to increase the characteristic electric 
field strength: the mode volume V may be decreased, the impedance n/e may be 
increased, or the group index n g may be brought close to zero. In the ideal white 
light cavity, n g — > 0, which implies that \£\ — > oo. 

One other interesting aspect of the ideal white light cavity is the mode energy. 
As n g — > from the positive side, U as defined in Eq. (6.12) also approaches zero. 
This is a result of a cancellation of energy between a negative stored energy (U s ) and 
a positive field energy (Uf). By analogy with the definition of Uf, we can define a 
field Hamiltonian (Hf) as 

H f = huj— fo f oH- l ) , (6.17) 

n g V 2J 

which gives a vacuum field noise energy of 

hion 
U. = ^-, (6.18) 
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R.V 



R<1 
W 



R=l 



Figure 6.1: A small cavity of length £ is coupled to a larger cavity of length L by 
a mirror with reflectance R. The small cavity contains a lossless, dispersive medium 
with a group index of n 9 . 



which becomes infinite for the closed cavity model as n g approaches zero. In atomic 
vapor systems currently used to achieve anomalous dispersion, n, e r , and \i r are very 
close to 1. The total vacuum field noise then reduces to 



U„ 



2n n 



(6.19) 



In the next section we will present an open cavity model of a white light cavity 
that will reveal the effect of the group index on the cavity noise spectrum and allow 
us to make specific calculations based on reasonable models for a white light cavity 
medium. 



6.3 Coupled cavity approach 

In the previous section, we saw that the standard expression for the electromagnetic 
energy density in a dispersive medium implies that dispersion may scale the ampli- 
tudes of the electromagnetic fields associated with frequency eigenmodes in a closed, 
lossless cavity. Because the closed cavity treatment neglects loss, the field modes 
in the model are delta functions in frequency and give no insight into the effect of 
dispersion on the power spectrum of the modes. For the case of a white light cavity, 
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the closed cavity treatment also suffers from the fact that the classical expression 
for the electromagnetic energy density of a white light mode suggests a total energy 
of zero. In this section we extend and corroborate the closed cavity approach by 
examining the effect of dispersive on the electromagnetic fields of a lossy cavity. 

Cavity loss means that the electromagnetic energy inside a cavity is coupled 
to some external system or systems. In practice, cavity modes are coupled both to 
external electromagnetic field modes and also to atomic systems which in turn may be 
coupled to other modes and systems. Electromagnetic field modes are equivalent to 
simple harmonic oscillators. In linear regimes, atomic systems may also be modeled 
as simple harmonic oscillators. Thus a simple and reasonable way to model loss is 
to couple the system explicitly to a set of simple harmonic oscillators. 

The simplest way to do this may be to couple the lossy cavity to a larger cavity. 
The lossy cavity mode then behaves as a coupled harmonic oscillator and the mode 
is referred to as a pseudomode. This lossy cavity pseudomode can be probed using 
the modes of the overall system, which behave as uncoupled harmonic oscillators, 
and are often called "true modes" or "universe modes [16]." A nice discussion of 
the interrelations between pseudomodes, quasimodes, and true modes may be found 
in a thorough paper by Dalton and colleagues [ ]. Our contribution here is to 
show not only that the true mode approach yields insight on the pseudomodes of a 
dispersive cavity, but that it does so in a way that not depend on the exact form of 
the electromagnetic energy density. 

Figure 6.1 depicts a simple 1 dimensional implementation of this approach in- 
volving three mirrors. The two outer mirrors are perfect, and the combined system 
is closed. The inner mirror has a reflectance R, is lossless, and is frequency indepen- 
dent for the range of frequencies that will interest us. These three mirrors form three 
distinct resonant systems: a longer reservoir cavity of length L, a shorter cavity of 
length £, and a closed system comprised of the these two cavities coupled together. 
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From here on we will refer to these resonant systems as the f-cavity, the L-cavity, 
and the c-cavity. The L-cavity is evacuated while the ^-cavity contains a dispersive 
but lossless medium with refractive index, relative permeability and relative permit- 
tivity all equal to 1 at line center. 5 Our object is to see how dispersion affects the 
pseudomodes of the £-cavity by using the modes of the c-cavity as probes. 

We will begin by using this model to corroborate the results of Section 6.2. By 
neglecting absorption/gain and group velocity dispersion, we will use the open cavity 
model rederive Eq. (6.19) without relying on any expression for the electromagnetic 
energy density in a dispersive medium. We will then extend this result by includ- 
ing group velocity dispersion using a specific model for an anomalously dispersive 
medium. 

6.3.1 Total field noise of an open dispersive cavity mode in 
the absence of loss, gain, and group velocity dispersion 

We we begin by using the model depicted in Figure 6.1 to calculate the total field 
energy in a dispersive cavity pseudo-mode. The total field energy (JJfi) of an £- 
cavity pseudo-mode is a combination of contributions from many c-cavity modes. 
Thus 

U ft = £AJ2 u f^ (6- 2 °) 

4>l = — 1T 

where A is the mode area, Ufa is the average energy density of the i'th c-cavity mode 
in the £-cavity, and the sum is taken across one £-cavity free spectral range. 

For the z'th c-cavity mode, the electric field amplitude in the ^-cavity (E#) is 



5 We choose n = 1 not only because it is close to the refractive index value associated 
with many possible white light cavity implementations, but also because it allows us to 
make a simple distinction between total and field electromagnetic energies. 
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related to the electric field in the L-cavity (E^) by 

E« = E Ll - J=—-, (6.21) 

where t is the transmission coefficient of the central mirror and <$& is the round trip 
phase of the ^-cavity at Wj. 

Using Eq. (6.21) and noting that the ratio of field energy densities is proportional 
to the squared modulus of the squared electromagnetic fields, we obtain 

Ufa _ 1 - R 



u L i l + J R-2 v A Rcos(0 i )' 
where we use the losslessness of the mirror to get tt* = 1 — R. 



(6.22) 



* =7r 1 — i? 

Un = IA Y" u Li = . (6.23) 

^ l + R-2VRcos(<j )l ) 

We know that if R were 1, the modes of the L-cavity would act as ordinary simple 
harmonic oscillators because we understand the quantization of electromagnetic fields 
in a vacuum. The c-cavity modes are perturbed by the effect of the small cavity. 
However, assuming that L is large and that this perturbation would therefore not 
change the simple harmonic oscillator character of the overall modes, we take the 
total vacuum energy of the overall cavity modes to be fru/2. 

To get an expression for uu, we divide by a weighted volume. If we take the mode 
area in each cavity to be A then the volumes of the L-cavity and the f-cavity are AL 
and AL The squared field in the f-cavity differs from the that in the L-cavityby 
the factor given in Eq. (6.22). The energy density in the small cavity differs from 
that of the field energy by the factor n g . Taken together, these considerations give 
an L-cavity energy density of 



hui l 

~2~LA + EA'' 



«« = V-TTT-STT. ( 6 - 24 ) 
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where 

£ = £n„ ^E . (6.25) 

Thus, the total field energy of the pseudo-mode may be written as 

Uft = l*Y —^^- a ^# • (6-26) 

As L grows large, the free spectral range of the c-cavity (Aui+e) falls below the 
minimum relevant resolution of the system and we may replace the sum over modes 
by a spectral integral without sacrificing accuracy Thus, 

' ~- w hw l 



U fe ^£A I — — ^^x 



l-R I A (6 ' 27) 

-dco. 



1 + R- 2VRcos(4>(co)) Aw l , 

We can calculate Aui+e by noting that the round-trip phase <fi for any cavity 
must change by 2tt between two adjacent resonances, giving 

A dd> Au 2 d 2 d> 
A0 = 2vr = n+1 - n = Aw/ + — - y^ + ... • 6.28 

du 2 du 2 

Since Aui+e is small for large L, Eq. (6.28) reduces to 

Auj L+i « 2iTdLo L+£ /d(f) L+e (6.29) 

when applied to the c-cavity. Using the definition of the group index and the c- 
cavity geometry, we find 

dfL+t 2L 2£ l-R 

= 1 n„ -= . (6.30) 

du L+l c c y l + R-2^Rcos((j) e ) 

Combining Eqs. (6.29) and (6.30) and using the definition 6.25 gives 

^-'"w^wr (6 ' 31) 
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Substituting Eq. (6.31) into Eq. (6.27) gives 

£ f* =7T fkJ 1 - R , n , 

U first— = du. 6.32 

1 7icJ^ = _ n 2 l+R-2y/RcOs(<f>(u)) 

We notice here that the term L + £'(u) resulting from the expression for dd> c /doj c 
cancels exactly with a similar effective length in the expression for the energy density. 
This is a particular example of what may be a more general correspondence between 
the derivative dd)/du and the electromagnetic storage capacity of lossless dispersive 
elements. 

To simplify Eq. (6.32) we now make the transformation 

du ->■ -£— dd> (6.33) 

to obtain 

U fi « * [^ - ^ #■ (6-34) 

To further simplify Eq. (6.34)) we make a third approximation by pulling oj and 
n 9 out of the integral. Doing so gives us 

tt hco ° l [^ 1 ~ R (a M 

Ufi rj / = dd), (6.35) 

n 2n g 27rJ (f)= _ ir l + R-2VRcos((j ) ) V ; 

and since 

1-R 



2tt, (6.36) 



- ff 1 + i?- 2^003(0) 

we find that 

^^g, (6.37) 

which is just a restatement of Eq. (6.19). 

The preceding derivation leading to Eq. (6.37) may be altered as follows to avoid 
reliance on any particular expression for the total electromagnetic energy density in 
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a dispersive medium. Equation (6.24) gives the energy density of a c-cavity mode 
in the L-cavityassuming a particular form for the energy density in the ^-cavity. 
Rewriting Eq. (6.24) without specifying a dispersive energy density gives 

hjJi X (a ?q\ 

Un = ; in—i (6.38) 

2 LA + lArf V ; 

where we make no assumptions on n, other than that it is finite. As L grows large, 
any finite value for n eventually becomes irrelevant: 

hm u Li = — — — . (6.39) 

L-s>oo Z LA 

Through an identical process, the free spectral range of the c-cavity, given by Eq. (6.29), 
becomes 

lim Au L+e = -^. (6.40) 

Limits (6.39) and (6.40) may be used in place of Eqs. (6.24) and (6.31) with no 
change to the overall argument, yielding the final result given by Eq. (6.37) without 
invoking a particular form the energy density of an anomalously dispersive medium. 

Eq. (6.37) is valid only where loss/gain and group velocity dispersion are both 
negligible over the cavity resonance. For media with substantial dispersion, this lim- 
its its applicability to narrow spectral ranges. For cavities exhibiting strong normal 
dispersion, this presents little difficulty because dispersion may narrow the resonance 
bandwidth of an already high-Q cavity. However, for the white light cavity in par- 
ticular, n 9 — > and the vacuum field energy predicted in Eq. (6.37) diverges. The 
open cavity model shows that this divergence comes is not due to higher buildup, 
but to a diverging bandwidth. In practice, the bandwidth of a given white light 
resonance is limited by group velocity dispersion (GVD) and by the cavity finesse. 
In the next subsection we explore the relationship between these two quantities in a 
specific white-light cavity realization. 
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6.3.2 Cavity noise in a symmetric gain-doublet system 

Causality dictates that every anomalously dispersive cavity will exhibit both group 
velocity dispersion (GVD) and loss and/or gain. These quantities could vanish at line 
center and so might be ignored if the resonance linewidth were sufficiently narrow. 
As an anomalously dispersive cavity more closely approaches a white light cavity, 
the cavity linewidth grows. A white light cavity would correspond to an infinitely 
wide linewidth. Any physically reasonable model of a white light cavity resonance 
must include GVD. 

A proper consideration of the total cavity noise must also account for the gain of 
the white light medium. Causality dictates that a medium that is both anomalously 
dispersive and transparent at a frequency u must exhibit gain somewhere not too 
distant from u . This gain may be difficult to maintain in a high-finesse cavity 
without inducing lasing. The difficulty may not be circumvented by making mirror 
reflectivity spectrally dependent without altering the anomalously dispersive nature 
of the cavity mode. This gain adds noise to our system. However, in the following 
calculations we will ignore its effect to show that anomalous dispersion alone is 
sufficient to require significant noise increases in a white light cavity. Because we 
will ignore gain in the paragraphs to come, the analysis contained in them will 
underestimate the total noise associated with the cavity mode. 

We will now solve for the pseudomode field energy given by Eq. (6.27) for a 
cavity that is homogeneously filled with a medium whose energy level diagram is 
depicted in Fig. 6.3. This medium produces a gain-doublet, as originally proposed 
by Steinberg and Chiao [24] , accomplished through a bichromatic Raman system as 
experimentally realized by Wang, Kuzmich, and Dogariu [8] . We take the two Raman 
gain lines to be centered about the frequency uj and separated by a distance of 2 A. 
The value for A may be tuned dynamically according to experimental expedience. 
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Figure 6.2: Energy level diagram for the Dichromatic Raman system. Pump lasers 
represented by lines Q\ and f2 2 provide symmetric gain lines around u . Detuning 
from | a) is large (A >> 2A). 



The linear susceptibility (x( w )) ma Y then be represented as [89, 83] 

XM = t , .-p + , a i -r ' 

u — UJ — A + zl LO — UJ + A + IL 



(6.41) 




-2 2 

w-u„ [MHz/2tt] 



Figure 6.3: The susceptibility of a typical bichromatic Raman system. The asym- 
metric line represents real susceptibility and the solid symmetric line represents the 
gain of the Raman system. The dashed line represents the total gain in the presence 
of a broadband absorber whose strength is chosen to match the Raman gain at line 
center. 
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9 10 11 12 
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Figure 6.4: Uf v is calculated for cavity lengths of lm (squares) and 0.25m (triangles) 
for varying values of finesse using Eq. (6.32) 



where T is the Raman transition line width, and the rates Mi and M 2 are given by 



Mi 



N(\u. ab ■ e\ 2 /2he )(\tt j \ 2 /A 2 ) } (j = 1,2) [ ]. These rates M x and M 2 depend 



on N, the number density of participating atoms, on fi a b ■ e, the dipole interaction 
between the exciting fields and the primary Raman transition, and on the Rabi 
frequencies of the two exciting fields. This last dependence means that they may 
be dynamically controlled by choosing the intensity of the pump lasers. Symmetry 
minimizes gain and group velocity dispersion at line center when the values of M% 
and M 2 are chosen such that M x = M 2 = M. For a given set of values u , A, and T, 
M may be chosen such that n g = at u . 



Figure 6.2 plots the susceptibility as a function of frequency for a typical bichro- 
matic Raman system where Mi = M 2 = M is chosen to minimize n g at line center. 
The dashed line indicates that this system may be combined with a broad-band ab- 
sorber to roughly eliminate the net gain at line center. This effect was approximated 
in the original experiment by Wang, Kuzmich, and Dogariu [ ]. Loss similarly played 
an important role in the more recent experiments reported by Pati et al. [31]. 
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Figure 6.4 shows total cavity vacuum field noise for increasing finesse values as 
calculated from Eq. (6.32). For both series of points, we used the values: T = 
1 x 10 6 Hz/2tt, A = 2 x 10 6 Hz/2tt, M « 3.45 x l(r 3 Hz/27r, u « 2.42 x 10 1 5Hz/27r, 
which were taken from the literature [31, 83] and then adapted to avoid negative 
group velocities near the gain lines (Negative group velocities near the gain lines 
lead to extra resonances in the wings. Causality dictates that these resonances exist 
at frequencies where the gain is higher or the loss lower than at line center and leads 
to unnecessarily large noise terms. In many cases, it is also expected to lead to 
lasing. This effect rules out the steady state stability of a large class of anomalously 
dispersive cavities.) We see that for both cavities the noise increases as higher finesse 
confines resonance to frequencies where group velocity dispersion has less effect. As 
finesse becomes infinite, so does Uf v despite the fact that Eq. (6.32) neglects any 
contribution of the gain of the active medium to the survival factor. 

6.3.3 Anomalous dispersion, spontaneous emission, and laser 
linewidth 

For a dipole emitter coupled to a large number of electromagnetic modes, the rate 
of spontaneous emission into one particular mode, for example the cavity mode, is 
proportional to the square of the electric field associated with that cavity mode. 
Then Eq. (6.15) suggests that the spontaneous emission rate into a particular cavity 
mode will scale with the squared electric field, which scales inversely with the group 
index. When a narrow-band approximation applies, this is the case. Thus, Eq. (6.37) 
provides a simple explanation for line narrowing Agarwal predicted in the case of 
lasing without inversion [ ]; the line width is already sufficiently narrow compared 
to the linewidth of the emitting particle such that further narrowing provides a simple 
effect on the spontaneous emission rate and the laser linewidth will scale with the 
group velocity. This result, although not explicitly applied in Agarwal's paper, can 
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explain his results. However, the result is not new. For example, the role of the 
group velocity in determining laser linewidth was understood by Henry [91], who 
framed the laser linewidth in terms of a ratio between the spontaneous emission rate 
and the number of photons in a mode and got an expression for the laser linewidth 
which is implicitly proportional to the group velocity. 

This same result applies to anomalously dispersive resonances, provided that the 
cavity linewidth is still narrow compared to the line shape of the gain medium. If 
this condition is not met, then the broadening due to anomalous dispersion is limited 
by the line shape of the gain medium. In any case, anomalous dispersion leads to an 
increase in laser linewidth just as it does in the linewidth of an empty cavity. 

The difference between our prediction for the effect of anomalous dispersion on 
laser linewidth and that of Shahriar et al. can be remedied if we redo their calculation 
while taking into account the effect of anomalous dispersion on the cavity photon 
lifetime. Following previous authors [92], they used a formula for the quantum limited 
linewidth that is equivalent to 

Aw, aser = — y/n, (6.42) 

where r c was a photon lifetime and n. They took r c to be the photon lifetime of the 
evacuated cavity. However, the photon lifetime for the anomalously dispersive cavity 
differs from that of the evacuated cavity in their case by the factor l/n 9 , where n 9 is 
a cavity averaged group index. If this extra factor is taken into account, then their 
prediction agrees with ours. 



6.4 Conclusion 

We have shown that the benefits achieved through intracavity anomalous dispersion, 
i.e. high bandwidth combined with high buildup, come only at the cost of increased 
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electromagnetic field noise. An infinite finesse combined with a perfect white light 
mode would lead to an infinite amount of field energy in the ground state of the cavity 
mode. Relaxing the infinite finesse condition to allow for cavity loss broadens the 
mode and relaxes the divergence, but still suggests a substantial increase in quantum 
field noise. 
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Summary and Outlook 



Give me a place to stand on, and I will move the Earth. 

-Archimedes 

According to Pappus [ ], Archimedes used the phrase in the above quote to 
express his enthusiasm for the principle of leverage. This dissertation is premised on 
the idea that the availability of controllable dispersion combined with only moderate 
absorption or gain provides us with a new physical lever. The range of physically 
viable group index values in relatively lossless media is now so extreme that it is 
hard to express this range succinctly. The group index can be made many orders 
of magnitude larger or smaller than its vacuum value. Its sign may be changed and 
many more orders of magnitude are accessible on the negative side. The question that 
we have posed is whether this lever can have interesting theoretical and experimental 
consequences. Chapters 4, 5, and 6 can be seen as explorations of this question. 
To answer these questions well, we had to come to a basic understanding of the 
meaning of dispersion. Hopefully, such an understanding was partially developed 
through Chapters 2 and 3. 
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Dispersion is, at its heart, a geometrical concept. The definition of the group 
velocity, v g = du/dk, is completely geometrical, and has geometrical consequences. 
In Chapter 2, we found that the spectrally averaged group index may be used to 
find the velocity of the planes of constant phase difference associated with copropa- 
gating plane waves. Since the position of a pulse is the consequence of interference 
between plane waves, the group index drives the velocity of the pulse. We joined with 
Rayleigh, Craven, and Candler of old in noting the important role that dispersion 
plays in spatial interferometers, and noted that dispersion alters the map between 
changes in frequency and changes in wavelength. We briefly overviewed some of the 
most basic techniques by which dispersion may be controlled and combined with 
limited loss and gain and vanishing group velocity dispersion for a narrow spectrum. 

Although it is possible to combine strong dispersion with minimal loss or gain, 
this is only possible to do for a narrow band of frequencies. For such a narrow 
band of frequencies and in the absence of permanent loss or gain, we found ourselves 
able to make simple statements about the relationships between dispersion and the 
distribution of energy between field and medium. In Chapter 3, we found that when 
strong dispersion is combined with negligible loss, the form of the electromagnetic 
energy density must be altered to account for dispersive energy exchange between 
the field and the medium. In the case of anomalous dispersion, the total energy 
is negative. We noted that this negative energy density is just a representation of 
energy that has left the medium but that will later return to it. Following Peatross, 
Ware, and Glascow [ ], we invoked the instantaneous spectrum as an explanation 
for this exchange. 

Energy exchange must be accompanied by momentum exchange. In Chapter 4 we 
saw that because of the effect of dispersion on the division of wave energy between 
field and medium, the nondispersive forms of the momenta associated with light 
propagating through a medium need to be adapted if dispersion is taken into account. 



110 



Chapter 7. Summary and Outlook 

We saw how the non-scattering force of an electromagnetic wave in a dispersive 
medium on an embedded particle may be partitioned into two components, one that 
is proportional to a temporal gradient and another that is proportional to a spatial 
gradient. By modulating the amount of dispersion in a system, we may modulate 
the ratio between temporal and spatial gradients, and so modulate the force and 
displacement of an embedded particle. 

Not only is it possible to set the spectral derivative of the refractive index to 
a particular value, but it is possible to alter this value in real time. In Chapter 5 
we explored the transformations wrought upon a pulse by a nonstationary medium 
that is also dispersive. We introduced a set of rules based on a simple symmetry 
that can be seen as a generalization of Snell's Law and of the reflecting Doppler 
effect. Using this set of rules, we showed that dispersion may be used to modu- 
late the frequency response of a pulse to a time dependent refractive index and the 
Doppler shift of a pulse reflecting off of a moving interface. We also showed how 
time-dependent dispersion may be coupled with spatially dependent dispersion to 
scale a given pulse longitudinally and temporally (and so spectrally) over several 
orders of magnitude without otherwise altering the shape of the pulse envelope. We 
then introduced specific boundary conditions that might be associated with slow 
changes to a medium and compared, for those boundary conditions, the effects of 
altering phase velocity in a nondispersive medium with the effects of altering group 
velocity in a dispersive medium on basic quantities like energy density, photon num- 
ber, Canonical momentum density and Abraham momentum density and listed the 
comparison in Table 5.1. 

In Chapter 6, we explored the effect of dispersion on the total vacuum field 
energy associated with a single cavity mode. We argued first on the basis of the 
dispersive energy density that was discussed in Chapter 3, and then based only on 
the geometric properties of a dispersive medium that the vacuum field energy is 
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modulated by dispersion using a modes-of-the-Universe approach. When the group 
index is increased by a factor of 2 and other optical quantities are left unchanged, 
the spectral width of the resonance decreases by a factor of 2 while the peak of 
the resonance is unchanged. The result is a 2-fold decrease in the total quantum 
field noise associated with that mode. In Appendix B, we extend this idea to other 
features of a cavity resonance. We see that in many respects the effect of scaling the 
group index associated with a particular cavity mode is exactly the effect that we 
would expect if we could scale the cavity length while preserving modal resonance 
at a particular frequency. 

As was mentioned in Chapter 1, Professor Robert Boyd's group at the University 
of Rochester is currently planning experiments to test some aspects of the theory 
developed in Chapter 4. The consideration of dispersion adds clarity to the ongoing 
discussion of optical momentum in a medium and I believe that there is room for 
additional contributions relating to the interplay between dispersion and momentum 
in a medium. When dispersion is considered, two experimentally relevant formula- 
tions for the momentum of a single photon appear, one that depends on the group 
index and the other that depends on the phase index. These two formulations for 
the momentum are related to two separable concepts, those of impulse transfer, and 
energy transport. 

It may be interesting to use dispersion as a lever to increase or dampen the 
frequency response of a pulse to the temporal changes in the refractive index of 
a medium as was proposed in Chapter 5. Because of the bandwidth limitations 
associated with non-lossy media exhibiting controlled dispersion, it is unlikely that 
temporal changes in the group velocity of a medium will lead to the generation of 
either extremely long or extremely short pulses. However, it should be feasible to 
change the bandwidth and longitudinal extent of a pulse over one or two orders of 
magnitude while leaving the pulse envelope otherwise unchanged. This would be 
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interesting to pursue experimentally. In addition, many of the results in Chapter 5 
are predicated on plane waves interacting with a surface that is sufficiently large 
and flat such that Snell's Law would hold if the surface were stationary For surfaces 
approaching the wavelength of light, it will be necessary to consider diffractive effects. 
Similar "temporally diffractive" effects are to be expected for interactions between 
surfaces and pulses that are limited in temporal extent. The spectral broadening 
associated with such effects may depend on the relative velocities of a given surface 
and wave in interesting ways. 

Finally, the effect of the dispersion on the quantum field noise associated with 
a single cavity mode needs to be checked experimentally, particularly for cavities 
where the average group index approaches zero. For this case, we predict enhanced 
spontaneous emission into the cavity mode so long as the empty cavity resonance 
is narrower than the free space line shape of the emitting particle. This could have 
practical application where there is a need for fast emission into a particular spatial 
mode and merits exploration. 
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Appendix A 



Forms of the group index 



Throughout this dissertation, many formulations of the n g are used. This appendix 
is dedicated to justifying and collecting these formulations. 

Although group velocity is not well defined over a finite spectrum in dispersive 
media, it has a sensible definition over an infinitesimal one, 

Thus, the lack of a definition over a finite spectrum is just a manifestation of the 
fact that v g is a function of frequency. The group index is defined, in analogy with 
the phase index, as [93] 

(A.2) 
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Expanding, we get 

dn , . rN 

n a = n + co-—. (A. 5) 

du 

We can rewrite the second term as a partial logarithmic derivative, giving 

dn 

n * = n+ dMu)- (A ' 6) 

One interesting thing about the logarithmic derivative is that constant multipliers are 
irrelevant. Thus, dln(u) = d\n(v) = dln(E) = drn(<7o), where E here refers to the 
energy per photon, Hlo. Another interesting thing is that an inversion corresponds 
to a sign change. Thus, dln(u) = — dln(Ao). 

This explains a few other ways of writing n g : 

(A.7) 

(A.8) 

(A.9) 

(A.10) 

Each of these relates simply to a derivative expression: 

n 9 = l~( nu ) ( A - n ) 

dv 

n g = ^(nE), (A. 12) 

n g = —— (na ), (A. 13) 

U(To 

n g = -Xl^-(n/X ). (A.14) 

(A. 15) 



All of these expressions are essentially in terms of the refractive index versus the 
photon (or polariton) energy. It is also illuminating to understand n g in terms of 
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what it does to the interaction between the photon energy and the wavelength. We 
begin by expanding dX/dv. 



dX 

dv 


d / c N 
dv \nv i 


Recalling 


Eq. A.ll, 


d\ 
dv 


cn g 
(nv) 2 


We rewrite this as 


n g = 


din X 
dlnv 



(A.16) 

This highlights the way in which dispersion alters the mapping between changes in 
spectral and spatial frequencies. 

The total phase associated with the propagation of a wave with spatial frequency 
k through a length £ of a homogeneous medium is given by <fi = kl. The total 
phase difference accrued by two neighboring frequencies gives another way to view 
the group index. 

4> dk ln g 

duj dio c 

Thus, 

A higher group index means that phase differences between frequency components 
increase more quickly with space. 

The quantity n g /n is in some sense a better measure of dispersion than n g itself. 
Since n g = n + codn/du, we may write 

n„ d In n , . 

-*■ = 1 + -tt— ■ A.18 

n dmoj 
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In an absorptionless, isotropic medium, n = y/e r u- r . Because a logarithmic deriva- 
tive does not change with its arguements are scaled by constants, we may then write 

n a d\n(en) 

— = 1+ ,, • 

n am u 

This may be expanded as 

n a ldlneldlnu , . „, 

9_ _1 j I <^_ l\ 1 Q\ 

n 2 dhiuj 2 dhi jj, 



Hi 



Appendix B 

Simulating a long cavity with a 
short one 



The optical characteristics of a cavity can be modified over a narrow frequency range 
upon the introduction of a medium that exhibits electromagnetically induced trans- 
parency (EIT). Cavity characteristics of free spectral range, resonance bandwidth, 
quality, photon lifetime, and mode volume all change in a consistent fashion. In each 
case, the property is changed as though the cavity were increased in length by a 
factor of the group index n g . With group indices of ~ 10 7 possible by EIT media, a 
100 yum cavity so filled would exhibit the optical properties of a 1 km cavity. 



B.l Introduction 

Many important parameters of an optical cavity, including free spectral range, res- 
onance bandwidth, quality, photon lifetime, and mode volume [94, 95], depend on 
the cavity length. We point out that when a cavity is filled with a medium that 
exhibits electromagnetically induced transparency (EIT) [ ], the vacuum versions 
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of the expressions for these properties for frequencies near EIT resonance must be 
altered in a way that amounts to introducing a new effective length of the cavity. 
Consistently, the new length becomes the evacuated length multiplied by the group 
index of the medium. The cavity properties change accordingly. Using EIT, re- 
searchers have generated group indices greater than 10 7 [ ]. A 100 /xm cavity filled 
with such a medium would then assume the characteristics of a 1 km cavity. In this 
Appendix, we introduce the relevant features of EIT and then discuss each of the 
mentioned parameters in turn, showing how the effective length in each case is the 
evacuated length scaled by the group index. 



B.2 EIT and the group index 

Optically, an ideal EIT 1 medium [ ] is like a dispersive vacuum. In a narrow fre- 
quency band about the EIT resonance, the absorption coefficient (a) and the re- 
fractive index (n) are near their vacuum values of and 1, respectively, while the 
group index (n 9 ) can be arbitrarily large [28, 7, 9]. The phase velocity approaches 
the vacuum speed of light (c) while the group velocity (v 9 ) can be much slower than 
c, according to v g = c/n g . 

A pulse flowing from a vacuum into an EIT medium will be spatially compressed 
along the axis of propagation by a factor n g . The envelope of the pulse will change 



X EIT is a phenomenon where a medium ordinarily opaque at a particular frequency can 
be made transparent on application of a laser at a second frequency. Typically, the effect is 
created by connecting three discrete levels through two one-photon transitions, each stimu- 
lated by a laser. A "probe" excites the normally opaque transition, and a "drive" completes 
an overall two-photon process to the third level. We consider the drive-beam frequency 
to be fixed and resonant with the transition it drives. Further, we take the medium as 
manipulated by the drive to be the EIT medium. Our concern is then with the properties 
of the EIT medium at the probe frequency. Resonance refers to the case where the probe is 
tuned to its corresponding one-photon resonance and the medium is transparent. Detuning 
will mean detuning of the probe beam only, and imperfect transparency. 
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abruptly in slope at the interface but will be otherwise continuous because of the 
identical refractive indices on either side of the boundary. No energy is reflected at 
the interface, for example; all is transmitted. That the pulse can be foreshortened 
and unchanged in field intensity is owed to the fact that an EIT medium stores 
and exchanges energy with the propagating electromagnetic field. The total energy 
density, /i, in an EIT medium is apportioned in general by [18, 10] 

H = Vm + Hf = n g fi f , (B.l) 

where [if is the energy density of the electromagnetic field, and [i m is the density 
of energy stored by the medium. Thus, an EIT-filled cavity may hold much more 
energy than is evident in the internal intensity alone. 

The foreshortening of the pulse highlights another effect of n g that is relevant 
to cavities. Pulse propagation can be seen as the propagation of a set of phase 
relationships between its frequency constituents. On entering a high n g medium, the 
pulse compresses because these phase relationships change n g times more quickly for 
each increment of displacement. Thus an EIT-filled cavity not only holds n g times 
more energy for a given [if but also brings about changes in phase relationships 
between different frequencies as though the cavity were n g times as long. These two 
related effects are at the heart of the ability of an EIT cavity of length L to mimic 
an evacuated cavity of length n g L. 



B.3 Free spectral range 

The free spectral range, Az/ FSR , of a cavity is the frequency difference between suc- 
cessive longitudinal cavity resonances and is usually given by [5, 6] 

Az/fsr = ^, (B.2) 
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where L is the round trip length of the cavity and n is presumed constant. To take 
the EIT medium into account, we constitute Ai/fsr afresh. 

A cavity resonance occurs when the total phase <fi accrued by light propagating 
over L is an integer multiple of 27r, or when kL = q27c, where k = 2-n j\ m is the 
angular wavenumber, X m is the wavelength of the light in the medium, and q is an 
integer. The difference in k between adjacent resonances at q + 1 and q is 27r/L, 
which in terms of frequency v = kc/(2irn) becomes 

v q+1 n(v q+1 ) - v q n{y q ) = c/L. (B.3) 

A Taylor expansion of Eq. B.3 with 

Vq+l = V q + AZ/ FSR 

/ \ / \ a & n 

™K+i) = n{v q ) + Az/ FSR — + . . . 

leads to the first-order approximation 

Az/fsr = -^y (B.4) 

n 9 L 

since n 9 {v q ) = v q n' (z/ g )+n(z/ g ). When n 9 — >■ n, as is ordinarily the case in transparent 
media, the more generally applicable Eq. B.4 reduces to Eq. B.2. For an empty cavity 
n = 1 and u q+ i exceeds v q by Az/fsr = c/L as seen from Eqs. B.2 and B.3. Use 
of the resonance condition and a vacuum dispersion relation of dcf)/dv = 2tt(L/c) 
yields the same result. When filled with an EIT medium, <fi instead changes by 
d^/du = 27in g (L/c). The difference z/ g+ i — v q becomes c/(n g L), and Eq. B.4 is 
recovered, which shows Az/fsr to be that of an empty cavity of L stretched by n g , a 
factor equal to the increased dispersion of the EIT medium. 
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B.4 Resonance bandwidth 

The resonance bandwidth, Au r , of a cavity is given by 

Au r = ^* (B.5) 

where F is the finesse of the cavity [5, 6]. The finesse depends only on mirror 
reflectivities and loss within the medium and is unaffected by changes in L to the 
extent that losses remain constant. Assuming a fixed finesse, we see from Eqs. B.4 
and B.5 that insertion of an EIT medium leads to 

Av T = ^, (B.6) 

n g L 

that is, the equivalently sharp Av r of an evacuated cavity of length n g L. 

The reduction of Av r , and consequent improvement in spectral discrimination, 
of a cavity through EIT was predicted by Lukin et al. 96]. Wang et al. \ 
later demonstrated a narrowing of Av r and suggested potential uses of intracavity 
EIT in areas of high-resolution spectroscopy, laser-frequency stabilization, and the 
generation of non-classical light. 



B.5 Quality 

The cavity quality Q for a resonant mode can be expressed by 

Q = v/Av r , 
which through use of Eq. B.6 leads to 

Q = F n -f (B.7) 

in the case of a cavity filled with an EIT medium. 
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The quality may also be cast as the ratio between the steady total cavity energy 
U and the cavity energy lost per optical cycle according to 



Q = 2nv 



-dU/dt 



From Eq. B.l, U ~ ngl^f for a cavity of fixed volume. The loss rate of energy, 
however, scales with \ij in the medium and at the mirrors, giving dU/dt ~ (Xf. The 
quality then increases by n g because the energy stored by the EIT medium cannot 
be dissipated or transmitted by mechanisms of optical loss, which determine dU/dt. 



B.6 Photon lifetime 

The lifetime r of a photon 2 in a cavity is proportional to Q, which by Eq. B.7 leads 
tor ~ n g L. A separate consideration of r better clarifies its association with the 
equivalent length n g L. 

It suffices to define r by the round-trip time for a photon multiplied by the 
expected number of circulations. For an evacuated cavity we have r = (L/c)(l — s) _1 , 
where s is the probability the photon will survive one round trip, and (1 — s) _1 is 
the number of round trips associated with a 1/e photon lifetime in a high-finesse 
cavity [ ]. Insertion of an ideal EIT medium will not change s, but will increase 
the round-trip time to give r = (n g L/c)(l — s) _1 . This time is equivalent to the 
round-trip time of an evacuated cavity with a physical length of n g L. 



2 For a highly dispersive medium like EIT, where the energy is no longer only in an 
electromagnetic form, we might replace the term photon with quantum of energy. 



124 



Appendix B. Simulating a long cavity with a short one 

B.7 Mode volume 

In an evacuated cavity, the mode volume V is the volume available for a photon 
mode and is proportional to the cavity length. We may figure the mode volume 
from U, the total energy within the mode, and /i, the energy density at a particular 
point, from V = U/u. [94, 95]. The mode volume finds practical use in determining 
the average field amplitude A ex y/U/V of a mode given the energy U it contains. 
When the cavity is filled with an EIT medium, however, we must bear in mind from 
Eq. B.l that /i is divided between /i m and /i/, leading us to change V and A from 
their empty-cavity expressions to V — n g U/u. and A ex ^U/(n g V). It is natural to 
factor the mode volume into components of area and length. In doing so, we are 
led to define a new effective mode volume V e jf = n g V whose area is unchanged but 
whose length is scaled by the group index. 

We may view the expanded V and reduced A of intracavity EIT in terms of where 
the photon, or equivalent quantum of energy, resides. When that energy takes the 
form of an electromagnetic field, we consider the field's spatial extent to define the 
mode volume. However, for n 9 -times longer than present as a photon, the energy 
is stored within the medium. The optical-field amplitude available to affect other 
atoms, for example to bring about the decay of excited atoms that do not participate 
in EIT, must be based on ^JU /n g , or equivalently a mode volume extended by n g . 



B.8 Caveats 

While the effective length n g L controls the cavity characteristics discussed, there 
are other considered uses and properties of a cavity for which the length remains 
L. We draw attention to two of these. First, the effective length does not offer the 
extrinsic benefits of physical size. A long cavity has an enhanced sensitivity, through 
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detectable changes to its FSR or finesse, to the refractive index or absorption of a 
trace gas that fills its volume. This advantage is not conferred on a short cavity 
through EIT. Second, the number of nodes in the cavity is unchanged by the EIT 
medium. Although EIT increases n g it does not alter n enough to disturb the node 
spacing of A m /2 significantly. 

EIT introduces a severe bandwidth constraint for high n g operation, with n g ^> 1 
and constant only over a narrow range of frequency 5v about the EIT resonance. 
Also, n g and Sis are in opposition; high n g associates with narrow Sis. Some addi- 
tional frequency range can be gained by detuning the EIT drive laser and accordingly 
detuning the probe to recover the overall two-photon resonance, a process that essen- 
tially shifts the center frequency of Sis. Large tuning ranges, however, would likely be 
obtained only by using different energy-level configurations in the same or another 
medium. 

Finally, in an EIT-enhanced cavity, a stable effective length will require both 
mechanical stability and a steady drive-laser intensity. A positive aspect of the latter 
dependence, however, is that one can scan the effective length of an EIT cavity over 
many orders of magnitude without using any moving parts. 3 



B.9 Conclusion 

With respect to many optical parameters, an EIT cavity of length L is identical to 
a vacuum cavity of length n g L. Because n g can be made extremely large we can 
simulate the properties of very large cavities with minute ones. We can narrow the 
resonance bandwidth well below that set by the empty-cavity finesse, for example. 
As n g is not difficult to control dynamically, it should be possible to sweep the 



3 The drive-laser intensity and n„ are inversely related, with n„ > 1 for EIT. 
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effective length of a cavity quickly and easily over a large dynamic range. It is also 
worth considering media beyond those that exhibit EIT, where n g can approach zero 
and even assume negative values. In these cases, artificial cavities are created with 
vanishing to negative effective lengths, something impossible to do with a vacuum 
cavity by any manipulation. On a practical level, the relation Az/fsr ~ n~ 1 points 
out a simple cavity-based approach to accurate measurements of n g . 
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